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Abstract: We use the AdS/CFT correspondence to study the thermalization of a 
strongly coupled conformal field theory that is forced out of its vacuum by a source 
that couples to a marginal operator. The source is taken to be of small amplitude 
and finite duration, but is otherwise an arbitrary function of time. When the field 
theory lives on i?'^^^'^, the source sets up a translationally invariant wave in the 
dual gravitational description. This wave propagates radially inwards in AdSd+i 
space and collapses to form a black brane. Outside its horizon the bulk spacetime 
for this collapse process may systematically be constructed in an expansion in the 
amplitude of the source function, and takes the Vaidya form at leading order in the 
source amplitude. This solution is dual to a remarkably rapid and intriguingly scale 
dependent thermalization process in the field theory. When the field theory lives on a 
sphere the resultant wave either slowly scatters into a thermal gas (dual to a glueball 
type phase in the boundary theory) or rapidly collapses into a black hole (dual to a 
plasma type phase in the field theory) depending on the time scale and amplitude 
of the source function. The transition between these two behaviors is sharp and can 
be tuned to the Choptuik scaling solution in R^'^. 
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1. Introduction 



The AdS/CFT correspondence identifies asymptotically AdS gravitational dynamics 
with the master field evolution of 'large A' field theories. In particular, it relates the 
evolution of spacetimes with horizons to the non equilibrium statistical dynamics of 
the high temperature phase of the dual field theory. This connection has recently 
been studied in detail in a near equilibrium limit. It has been established that the 
spacetimes that locally (i.e. tube wise) approximate the black brane metric obey the 
equations of boundary fluid dynamics with gravitationally determined dissipative 
constants. ^ The equations of fluid dynamics are thus embedded in a long distance 
sector of asymptotically AdS gravity, a fascinating connection that promises to prove 
useful in many ways. 

Given the success in using gravitational physics to study near equilibrium field 
theory dynamics, it is natural to attempt to use gravitational dynamics to study far 
from equilibrium field theory processes. In this paper we will study the gravitational 
dual of the process of equilibration; i.e the dynamical passage of a system from a 



iSee § I, I i I i g I 10 111, 111, 111 

understanding of this connection. See the reviews pi 
work. See also M |25|, |2^, |27|, M |2g|, ISCl, ISll, |32l pi bI 



111, 111, for a recent structural 
0, 0for references to important earlier 
for related line of development. 
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pure state in its 'low temperature' phase to an approximately thermalized state in its 
high temperature phase (see [||, |6|, g^, |3|, |9|, H, |T], |4|, ||, ||, |5H|, g^, Hf for 



closely related earlier work and ^ ^ for analyses of thermalization directly 

in large N gauge theories). As has been remarked by several authors, this process is 
dual to the gravitational process of black hole formation via gravitational collapse. 
The dynamical process is fascinating in its own right, but gains additional interest in 
asymptotically AdS spaces because of its link to field theory equilibration dynamics. 
In this paper we study asymptotically AdS (and briefly asymptotically flat) collapse 
processes in a weak field limit that displays rich dynamics while allowing for analytic 
control. 

An AdS collapse process that could result in black hole formation may be set 
up, following Yaffe and Chesler [Q, as follows . Consider an asymptotically locally 
AdS spacetime, and let TZ denote a finite patch of the conformal boundary of this 
spacetime. We choose our spacetime to be exactly AdS outside the causal future 
of TZ. On TZ we turn on the non normalizable part of a massless bulk field. This 
boundary condition sets up an ingoing shell of the corresponding field that collapses 
in AdS space. Under appropriate conditions the subsequent dynamics can result in 
black hole formation. 

In this paper we will study the AdS collapse scenario (plus a fiat space counter- 
part) outlined in the previous paragraph in a weak field limit; i.e. we always choose 
the amplitude e of the non normalizable perturbation to be small. In the interest of 
simplicity we also focus on situations that preserve a great deal of symmetry, as we 
explain further below. In the rest of this introduction we describe the three classes 
of collapse situations we study, and the principal results of our analysis. ^ 



1.1 Translationally invariant asymptotically AdSd+i collapse 

In the first part of this paper we analyze spacetimes that asymptote to Poincare patch 
AdSd+i space and turn on non normalizable modes on the boundary ^. We choose 
our non normalizable data to depend on the boundary time but to be independent 
of boundary spatial coordinates. Moreover, our data has support only in the time 
interval v G (0, 6t), i.e. our forcing functions are turned on only over a limited time 
interval. Our boundary conditions create a translationally invariant wave of small 
amplitude e near the boundary of AdS. This wave then propagates into the bulk of 
AdS space. 



^In most of the bulk of the text of this paper we only present formulae for asymptotically AdSd+i 
spacetimes for the smallest nontrivial value of d namely d = 3. As we explain in the Appendix ^ 
however, most of the qualitative results of our analysis apply to arbitrary odd d for rf > 3 and also 
plausibly to arbitrary even d for d > 4. 

for other work on Poincare patch AdS solutions forced by time dependent 
non normalizable data 
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Figure 1: Cross section of the causal diagram for the collapse process in an asymptotically 
global AdSd+i space. The conventional Penrose diagram for this process would include only 
the half of the diagram that to the right of its vertical axis of symmetry 



In section ^ and Appendices ^ and [B.l| of this paper we demonstrate that this 



wave always results in black brane formation at small amplitude (see figure |I] for the 
Penrose diagram of the analogous process in an asymptotically global AdS space). 
Outside the event horizon, this black brane formation process is reliably described by 
a perturbation expansion in the amplitude. At leading order in perturbation theory 
the spacetime set up by this wave takes the Vaidya form^ (||59|, 0, see e.g.[|62 
for a review) 



M{v] 



ds' = 2drdv - [r' - ) dv' + r'dxi. (1.1) 



This form of the metric is exact for all r when f < 0, and is a good approximation 

2 

to the metric for r 3> ^-jj- when f > 0. Our perturbative procedure determines 
the function M{v) in ( p..l|) in terms of the non normalizable data at the boundary; 



4 



The Vaidya metric is an exact solution for the propagation of a null dust - a fluid whose stress 



tensor is proportional to pk^k^, for a lightlike vector fc^t (k^ = dr in (1.1)). Note that pk^^ki, is 
also the stress tensor of a massless field in the eikonol or geometric optics approximation. The 
AdS- Vaidya metric has been studied before in the context of the AdS /C FT correspondence in, for 
instance, [pSl. 
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M(f ) turns out to be of order ^. M{v) reduces to constant M for v > 6t in odd 
d and asymptotes to that value (like a power in ^) in even d ^ In either case the 
spacetime (|1 . 1| ) describes the process of formation of a black brane of temperature 



2 

T ~ ^ over the time scale of order 6t. Note that the time scale of formation of the 
brane is much smaller than its inverse temperature. This fact allows us to compute 
the event horizon of the spacetime (|1 . 1| ) in a simple and explicit fashion in a power 

2 

series in 6tT ~ e^. To leading order in e the event horizon manifold is given by 



rniv) = Md V > 

Ml (1-2) 

rH{v) = — V < 

1 - Mdv 



All of the spacetime outside the event horizon ( [L.2D lies within the domain of validity 
of our perturbative procedure. Of course perturbation theory does not accurately 
describe the process of singularity formation of the black brane. However the region 
where perturbation theory breaks down (and so (|1 . 1|) is not reliable) is contained 
entirely within the event horizon of (|1 . 1| ) . Consequently, the region outside pertur- 
bative control is causally disconnected from physics outside the event horizon, so 
our perturbation procedure gives a fully reliable description of the dynamics outside 
the event horizon. It follows in particular that any singularities that develop in our 
solution are is always shielded by a regular event horizon, in agreement with the 
cosmic censorship conjecture. 

In section |^ we demonstrate that the corrections to the Vaidya metric (|1 . 1| ) may 
be systematically computed in a power series in positive fractional powers of e. At 
any order in the perturbation expansion, the metric may be determined analytically 
for times v <^ T^^ (T is the temperature of the eventually formed brane). At times 
of order or larger than T~^, perturbative corrections to the metric are determined in 
terms solutions of universal (i.e. independent of the form of the perturbation) linear 
differential equations which we have only been able to solve numerically. Even at late 
times, however our perturbative procedure analytically determines the dependence 
of observables on the functional form of the non normalizable perturbation, allowing 
us to draw conclusions that are valid for small amplitude perturbation of arbitrary 
form. 

Let us now word our results in dual field theoretic terms. Our gravity solu- 
tion describes a CFT initially in its vacuum state. Over the time period (0, St) the 



^More precisely, let (/>o(i') = e xijt) where x is a function that is defined on (0,1). Then the 
energy of the resultant black brane is ^^p- x A[x] where A[x] is a functional of x(a;) that is computed 
later in this paper. 

^M{v) is defined as the coefiicient of the term in the metric, in an expansion around small 
r. In even d this turns out not to be equal to the mass density of the system, the coefficient of the 
same term in the metric when expanded around large r. 

"^We thank M. Rangamani for discussions on this point. 
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field theory is perturbed by a translationally invariant time dependent source, of 
amplitude e, that couples to a marginal operator. This coupling pumps energy into 
this system. Our perturbative gravitational solution gives a detailed description of 
the subsequent equilibration process; in particular it gives a precise formula for the 
temperature of the final equilibrium configuration as a function of the perturbation 
function. It also, very surprisingly, asserts that for some purposes^ our system ap- 
pears to thermalize almost instanteneously at leading order in e. We pause to explain 
this in detail. ^ 

A field theorist presented with a fiow towards equilibrium might choose to probe 
this fiow by perturbing it with an infinitesimal source, localized at some time. He 
would then measure the subsequent change in the solution in response to this pertur- 
bation. However note that the spacetime in ( |1.1|) is identical to the spacetime outside 
a static uniform black brane for v > 6t when d is odd (and for v ^ 6t for even d). 
It follows that the response of our system to any boundary perturbation localized at 
times V > 6t in odd d (and aX v ^ 6t in even d) will be identical to the response of a 
thermally equilibrated system to the same perturbation. In other words our system 
responds to perturbations at v > 6t as if it had equilibrated instanteneously. 

A field theorist could also characterize a fiow towards equilibrium by recording 
the values of all observables as a function of time (in the absence of any further per- 
turbation). The full set of observables consists of expectation values of the arbitrary 
product of 'gauge invariant' operators, i.e. quantities that in a gauge theory would 
take the form 



{TrOiTr02...TrOn). 

In this paper we work in the strict large N limit (i.e. the strictly classical limit from 
the dual bulk viewpoint). In this limit trace factorization (or the classical nature 
of the dual bulk theory) ensures that the expectation value of products equals the 
product of expectation values. In other words our set of observables is given precisely 
by the one point functions of all gauge invariant operators. 

Now note that expectation values of all local boundary operators are determined 



by the bulk solution in a neighborhood of the boundary values. As the metric (|1.1| ) is 
identical to the metric of a uniform black brane in the neighbourhood of the boundary 
when V > 6t, it follows that the expectation value of all local boundary operators 



^In particular, in even bulk space time dimensions, one point functions of all local operators 
reduce to their thermal values as soon as the perturbation is switched off. While thermalization of 
one point functions is not instantaneous in odd bulk space times it appears to take place over a 
time scale of order St. 

^The rest of this subsection was worked out in collaboration with O. Aharony, B. Kol and S. 
Raju. See also the paper |Q, by Lin and Shuryak, for a very similar earlier discussion. We thank 
E. Shuryak for bringing this paper to our attention. 
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reduce instantaneously to their thermal values in odd d (and when v ^ 6t in even 
d). Consequently, all local operators appear to thermalize instanteneously. 

Not all gauge invariant operators are local, however. A field theorist could also 
record the values of non local observables, like circular Wilson Loops of radius a, 
as a function of time. As nonlocal observables probe the spacetime away from the 
boundary, their expectation values reduce to thermal results only after a larger time 
that depends on the size of the loop (this time is proportional to a at small a). 
So a diligent infinite field theorist would be able to distinguish (|1 . 1|) from abso- 
lute thermal equilibrium at times greater than 6t, but only by keeping track of the 
expectation values of non local observables. 

If one were to retreat away from the large A^ limit one would find large new classes 
of gauge invariant observables; the connected correlators of, for instance, local gauge 
invariant operators. Such correlators also sample spacetime away from the boundary, 
the distance scale of this nonlocal sampling being set by the separation between 
the operator insertions (see |]SS[ for a detailed discussion of properties of correlation 
functions in asymptotically AdS Vaidya type metrics). As in our discussion of Wilson 
loops above, the time scale for thermalization of such connected correlators is set by 
their separation (it is proportional to their separation when this separation is small). 

As we have seen, the time scale of equlibration of the solutions described in 
this paper depend on the precise question you ask about it. We would now like to 
describe a concrete and possibly practically important experimental sense in which 
our system behaves as if it were instanteneously thermally equilibrated. 

Consider the response of a CFT in its vacuum to a forcing function that varies 
- though only slowly - with x. We anticipate that at v = 6t the corresponding 
spacetime is locally (tube wise) well described by a black brane metric with a value 
of the temprature that varies with x (see ( ^.1|) and the discusson arounf it in section 
I). According to i |, |, i |, i 0, i |, H [lT|, O, ^ [l. Hi d, 0, H m, 10], 
the subsequent evolution of our system is governed by the equations of boundary 
fiuid dynamics. The initial conditions for the relevant fiuid fiow are given at f = 6t. 
Consequently an experimentalist who observes the subsequent fiuid fiow, and back 
calculates, would conclude that his system was thermahzed at v = 6t. 

The thought experiment of the previous paragraph is reminiscent of situation 
at the RHIC experiment. The back calculation described in this paragraph, in the 
context of that experiment, suggests that the RHIC system is governed by fiuid 
dynamics at times of order 0.5 fermi after the collision, much faster than suggested 
by naive estimates for thermalization time (see |^3|] and references therein). It is 
natural to wonder whether the mechanisims for rapid equilibration of this paper have 
qualitative applicability to the RHIC experiment. We leave a serious investigation 
of this question to future work. 

In summary, ( |1 . 1| ) describes a system whose response to additional external per- 
turbations at V > 6t is identical to that of a thermally equilibrated system and whose 
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one point functions of local operators also instanteneously thermalize. However ex- 
pectation values of non local observables (or correlators) thermalize more slowly, 
over a time scale that depends on the smearing size of the observable (or correlator) . 
We find instantenous thermalization of expectation values local operators and the 
scale dependence in the process of equihbration fascinating. In fact this discussion 
is reminiscent of precursors in the AdS / CFT correspondence [0, ^ . 

We emphasize that our discussion of thermalization applies only at leading order 
in e expansion. Indeed our analysis was based on ( [L . 1|) which accurately describes 
our spacetime only at leading order in e. At sub leading orders ( |1.1|) is corrected by 
perturbations that decay to the black brane result only over the time scale 1/T, in 
accordance with naive expextation. Consequently, the instantaneous thermalization 
of expectation values of local operators is corrected by sub leading equilibration pro- 
cess that take place over the time scale 1/T, the thermalization of linear fluctuations 
about a brane of temperature T. Note, in particular, that we have no reason so 
suspect that thermalization occurs over a time period that is faster than the naive 
estimate v = when e is of order unity or larger. 

1.2 Spherically symmetric collapse in fiat space 

We next turn to the perturbative study of spherically symmetric collapse in an asymp- 
totically fiat space. Consider a spherically symmetric shell, propagating inwards, fo- 
cused onto the origin of an asymptotically fiat space. Such a shell may qualitatively 
be characterized by its thickness and mass, or (more usefully for our purposes) by 
the Schwarzschild radius associated with this mass. It is a well appreciated fact 
that this collapse process may reliably be described in an amplitude expansion when 
y = ^ is very small. The starting point for this expansion is the propagation of a 
free scalar shell. This free motion receives weak scattering corrections at small y, 
which may be computed perturbatively. 

In section | of this paper we demonstrate that this fiat space collapse process 
may also be reliably described in an amplitude expansion at large y. In section ^ 
and Appendix p.2| we study this collapse process mainly in odd d (i.e. in even bulk 
spacetime dimensions). The starting point for this expansion is a Vaidya metric 
similar to ( |1.1| ), whose event horizon we are able to reliably compute in a power 
series expansion in inverse powers of y. Outside this event horizon the dilaton is 
everywhere small and the Vaidya metric receives only weak scattering corrections 
that it may systematically be computed in a power series in - at large y. As in the 
previous subsection, our perturbative procedure is not valid everywhere; however the 
breakdown of perturbation theory occurs entirely within the event horizon, and so 
does not impinge on our control of the solution outside the event horizon. 

At early times we are able to determine the perturbative corrections to the metric 
(order by order in ^) in an entirely analytic manner. However late time corrections 
to the metric are computed in terms of the solutions to relevant universal linear 
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differential equations, which we have not been able to solve analytically. However 
our perturbative solutions carry a considerable amount of information, even in the 
absence of an explicit analytic solution to the relevant differential equation. As an 
example, in section ^ we determine the fraction of energy of the incident pulse that 
is radiated back out to infinity to nontrivial leading order in the expansion in -. We 
are able to analytically determine the dependence of this fraction on the shape of the 
incident pulse upto an overall constant (see ( ^.34| )). The determination of the value of 
this constant requires knowledge of the explicit solution of the 'universal' differential 
equation listed in section ^, and may presumably be determined numerically. 

An order parameter (the presence of an event horizon at late times) distin- 
guishes small y from large y behavior, so the transition between them must be sharp. 
This observation was originally made about twenty years ago in classic paper by 
Christodoulou (see ||68l) and references therein) who rigorously demonstrated that 



collapse at arbitrarily large y results in black hole formation, while collapse at small 
y does not. As the fascinating transition between small and large y behaviors (which 
has been extensively in a programme of numerical relativity initiated by Choptuik 
^ ) ^°) presumably occurs at y of order unity. Consequently it cannot be studied 



in either the small y or the large y expansions described in our paper. 

As we do not have a holographic description of gravitational dynamics in an 
asymptotically fiat space, we are unable to give a direct dual field theoretic interpre- 
tation of our results reviewed in this subsection. See however, the next subsection. 

1.3 Spherically symmetric collapse in asymptotically global AdS 

The process of spherically symmetric collapse in an asymptotically global AdS space 
constitutes an interesting one parameter interpolation between the collapse processes 



described in subsections |LT] and [L^. We study such collapse processes in section ^ 
of this paper. In section ^ we have studied this collapse situation in detail only in 
d = 3. In this subsection we report the generalization of these results to arbitrary 
odd dimension, which may qualitatively be inferred from the results of Appendix |B[ 
Consider a global AdS space, whose boundary is taken to be a sphere of radius 
R X time. Consider a collapse process initiated by radially symmetric non normal- 
izable boundary conditions that are turned on, uniformly over the boundary sphere, 
over a time interval 6t. The amplitude e of this perturbation together with the di- 
mensionless ratio x = ^, constitute the two qualitatively important parameters of 
this perturbation. In the limit x — it is obvious that the collapse process of this 
subsection effectively reduces to the Poincare patch collapse process described in 
subsection p..l| , and results in the formation of a black hole that is large compared to 
the AdS radius (and so locally well approximates a fiat black brane); quantitatively 



i°See 1^, |7l| for reviews and |7|, ^ ^ for recent work interpreting this transition 

in the context of the AdS/CFT correspondence. 
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this turns out to work for a; -C e^. When x ^ ea the most interesting part of the 
collapse process takes place in a bubble of approximately flat space. In this case the 
solution closely resembles a wave propagating in AdS space at large r, glued onto 
a flat space collapse process described in subsection ^]2[ Following through the 
details of the gluing process, it turns out that the inverse of the effective fiat space y 

2d-2 

parameter (see subsection |1.2|) is given by ^ ^ . The parameter y is of order unity 



when X ~ e''-!. We conclude that the end point of the global AdS collapse process 
is a black hole for x -C e<*^ but a scattering dilaton wave for x ^ e^^. 

d-2 

The minimum mass of black holes formed through this process is (we work 
in units in which the mass of the black hole is simply the long time value of the 
parameter M in ( [4.10[ ), the global analogue of Let us contrast this with the 



minimum mass of black holes that we expect to be produced when we pump energy 
into the more slowly (i.e. through a forcing function whose time variation is of order 
-^) but over a long time period. As we have described above, slow forcing deposits 
energy into the gravitational thermal gas. By continually forcing the system one 
creates a thermal gas of increasing energy. At a critical energy density of order 
^, however, density fluctuations in this thermal gas become unstable |]78|; the end 



point of this instability is believed to be a black hole. Clearly this slow pumping 
in of energy produces black holes of energy or greater. It follows that black hole 
production can be produced more efficiently (i.e. at lower energies) via rapid forcing 
than via a slow pumping in of energy into the system. 

As we have explained above, when e ^ 1 and when x <^ e^^, we are able to 
reliably establish black hole formation within perturbation theory (see figure |1] for a 
Penrose diagram of this process). As in the previous two subsections, the starting 
point of the perturbative expansion always turns out to be a metric of the Vaidya 
form, whose event horizon we are able to reliably compute. Our metric receives only 
small scattering corrections outside the event horizon. Although the perturbative 
procedure breaks down badly near the black hole singularity, that is irrelevant for 
the physics outside the event horizon. 



-'^^This statement is only correct at times u <C To see why recall that when a collapsing shell 
in flat space forms a black hole, some of its energy is radiated out to I+. The resolution of the 
infalling shell into a static black hole and plus a shell radiated out to infinity occurs over a time 
scale set by rn the Schwarszshild radius associated with the infalling matter. In AdS space this 
shell eventually reflects off the boundary of AdS space at times of order (note that this is a much 
larger length scale than rn when the black hole is small enough) and then is refocussed on the 
origin of space. This process repeats itself unendingly; evcntauUy all of the energy of the intitial 
shell is absorbed by the black hole. Consequently AdS collapse processes always differ significantly 
from their flat space counterparts for v ^. In particular, while such a process can result in the 
formation of arbitrarily small mass black holes over time scale , the mass of black holes created at 
long times is bounded from below (see below for an estimate). We thank V. Hubeny for a discussion 
on this point. 
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On the other hand, when x 3> (but at small e), the incident waves simply 
scatter through the origin, and subsequently undergo periodic motion in AdS space. 
This free motion is corrected by interaction effects that will eventually cause this 
dilaton pulse to deviate significantly from its free motion over a time scale that we 

2ci-2 

expect to scale like a positive power of ^ 2 times the inverse radius of the sphere 

e d-2 

12 



Thermal 




Amplitude (epsilon) 

Figure 2: The 'Phase Diagram' for our dynamical stirring in global AdS. The final 

2 

outcome is a large black hole for x <C (below the dashed curve), a small black hole for 

1 1 
X <C e<*-i (between the solid and dashed curve) and a thermal gas for x 3> e<*-i . The solid 

curve represents non analytic behaviour (a phase transition) while the dashed curve is a 

crossover. 

Let us now reword our results in field theory terms. Any CFT that admits a two 
derivative gravity dual description undergoes a first order finite temperature phase 
transition when studied on S'^~^. The low temperature phase is a gas of 'glueballs' 
(dual to gravitons) while the high temperature phase is a strongly interacting, dis- 
sipative, 'plasma' (dual to the black hole). The gravitational solutions of this paper 
describe such a CFT on a sphere, initially in its vacuum state. We then excite the 
CFT over a time 6t by turning on a spherically symmetric source function that cou- 
ples to a marginal operator. The most important qualitative question about the 
subsequent equilibration process is: in which phase does the system eventually settle 
down within classical dynamics (i.e. ignoring tunneling effects) ? Our gravitational 

solutions predict that the system settles in its free particle phase when x e~ 

1 

but in the plasma phase when x ^ e^^-i. As in subsection p..l| the equilibration in 

the high temperature phase is almost instantaneous. However equilibration in the 

"'^^We expect this pulse to thermalize over an even longer time scale, one that scales as a positive 
power of the larger number, 



- 11 - 



low temperature phase appears to occur over a much longer time scale. We note 
also that the transition between these two end points appears to be singular (this 
is the Choptuik singularity) in the large limit. This singularity is presumably 
smoothed out by fluctuations at finite N, a phenomenon that should be dual to the 
smoothing out of a naked gravitational singularity by quantum gravity fiuctuations. 

In the rest of this paper and in the appendices we will present a detailed study of 
the collapse scenarios outlined in this introduction. In the last section of this paper 
we also present a discussion of our results. 

2. Translationally invariant collapse in AdS 

In this section we study asymptotically planar (Poincare patch) AdSd+i solutions to 
negative cosmological constant Einstein gravity interacting with a minimally coupled 
massless scalar field (note that this system obeys the null energy condition). We focus 
on solutions in which the boundary value of the scalar field takes a given functional 
form 0o(^) ill the interval {0,6t) but vanishes otherwise. The amplitude of (poiv) 
(which we schematically refer to as e below) will be taken to be small in most of this 
section. The boundary dual to our setup is a d dimensional conformal field theory on 
perturbed by a spatially homogeneous and isotropic source function, 4>o{v), 
multiplying a marginal scalar operator. 

Note that our boundary conditions preserve an R^'^ x SO{d — 1) symmetry (the 
W^^^ factor is boundary spatial translations while the SO{d— 1) is boundary spatial 
rotations). In this section we study solutions on which W^^^ x SO{d — 1) lifts to 
an isometry of the full bulk spacetime. In other words the spacetimes studied in 
this section preserve the maximal symmetry allowed by our boundary conditions. 
As a consequence all bulk fields in our problem are functions of only two variables; 
a radial coordinate r and an Eddington Finkelstein ingoing time coordinate v. The 
chief results of this section are as follows: 

• The boundary conditions described above result in black brane formation for 
an arbitrary (small amphtude) source functions 4>o{v). 

• Outside the event horizon of our spacetime, we find an explicit analytic form 
for the metric as a function of 0o(^)- Our metric is accurate at leading order 
in the e expansion, and takes the Vaidya form (|1 . 1| ) with a mass function that 
we determine explicitly ClS db function of time. 

-^^See section ^ for a discussion of the effects of potential Gregory-Laflamme type instabilities 
near this singular surface. 
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• In particular, we find that the energy density of the resultant black brane is 
given, to leading order, by 

^-A((S)7l((™>)i 

in odd d and by 

(ft If i+l d+l 

^^ = - ^d_i)2d J dhdhd,^- Mti)Hti-t2)e{h-h)d,^' Mt2) (2.2) 

2 

in even d. Note that, in each case, C2 ~ jf^yt- 

• We find an explicit expression for the event horizon of the resultant solutions, at 
leading order, and thereby demonstrate that singularities formed in the process 
of black brane formation are always shielded by a regular event horizon at small 
e. 

• Perturbation theory in the amplitude e yields systematic corrections to this 
leading order metric. We unravel the structure of this perturbation expansion 
in detail and compute the first corrections to the leading order result. 

While every two derivative theory of gravity that admits an AdS solutions admits 
a consistent truncation to Einstein gravity with a negative cosmological constant, the 
same statement is clearly not true of gravity coupled to a minimally coupled massless 
scalar field. It is consequently of considerable interest to note that results closely 
analogous to those described above also apply to the study of Einstein gravity with 
a negative cosmological constant. In Appendix ^ we analyze the process of black 
brane formation by gravitational wave collapse in the theory of pure gravity (similar 
to the set up of |S^), and find results that are qualitatively very similar to those 
reported in this section. The solutions of Appendix ^ yield the dual description of a 
class of thermalization processes in every 3 dimensional conformal field theory that 
admits a dual description as a two derivative theory of gravity. In fact, the close 
similarity of the results of Appendix ^ with those of this section, lead us to believe 
that the results reported in this section are qualitatively robust. In particular we 
think it is very likely that results of this section will qualitatively apply to the most 
general small amplitude translationally invariant collapse process in the systems we 
study. 

2.1 The set up 

Consider a minimally coupled massless scalar (the 'dilaton') interacting with negative 
cosmological constant Einstein gravity in d + 1 spacetime dimensions 

S = l d'^^^xVg (i? - - (2.3) 
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The equations of motion that follow from the Lagrangian ( |2.3|) are 



Ef,^ = G^^ - -d^(j)d^(f) + + -^{d(j)y ^ 

vV = o 

where the indices /i, u range over all d + 1 spacetime coordinates. As mentioned 
above, in this section we are interested in locally asymptotically AdSd+i solutions to 
these equations that preserve an R^'^ x SO{d— 1) symmetry group. This symmetry 
requirement forces the boundary metric to be Weyl flat (i.e. Weyl equivalent to flat 
^d-i.i-j. ]2owever it allows the boundary value of the scalar field to be an arbitrary 
function of boundary time v. We choose this function as 

Mv)=0 {v<0) 

M-v) < e {0<v<6t) (2.5) 
00 (t;) =0 {v> 6t) 

(we also require that (f)o{v) and its first few derivatives are everywhere continuous. ^^). 

Everywhere in this paper we adopt the 'Eddington Finkelstein' gauge Qrr = 
Qri = and Qrv = 1- In this gauge, and subject to our symmetry requirement, our 
spacetime takes the form 

ds'^ = 2drdv — g{r, v)dv'^ + /^(r, v)dx1 

(2.6) 

= 0(r, v). 

The mathematical problem we address in this subsection is to solve the equations 
of motion ( p.4| ) for the functions 0, / and subject to the pure AdS initial conditions 

g{r-, v) = (t> < 0) 

/(r,t.)=r (t;<0) (2.7) 
0(r,t;) = O (w < 0) 

and the large r boundary conditions 

hm ^ = 1 

hm ^ = 1 (2-8) 
lim 0(r, v) = 0o(^) 

r— >oo 

The Eddington Finkelstein gauge we adopt in this paper does not completely fix 



gauge redundancy (see [53| for a related observation). The coordinate redefinition 



"'^^We expect that all our main physical conclusions will continue to apply if we replace our (j)o - 
which is chosen to strictly vanish outside (0, 6t) - by any function that decays sufficiently rapidly 
outside this range. 
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r = f + h{v) respects both our gauge choice as well as our boundary conditions. In 
order to completely define the mathematical problem of this section, we must fix this 
ambiguity. We have assumed above that /(r, v) = r + 0{l) at large r. It follows that 
under the unfixed diffeomorphism, f{r,v) f{r,v) + h{y) + 0{l/r). Consequently 
we can fix this gauge redundancy by demanding that /(r, v) ^ r + 0{l/r) at large r. 
We make this choice in what follows. As we will see below, it then follows from the 
equations of motion that g{r) = + 0{1). Consequently, the boundary conditions 
]8D on the fields g, f and (p, may be restated in more detail as 

g{r,v) = r' (^1 + 

f{r,v) = r(^l + 0{^)^ (2.9) 
r 



Equations ( |2.4| ), ( p.6|) , ( |2.7| ) and ( p.9D together constitute a completely well defined 



dynamical system. Given a particular forcing function (poiv), these equations and 
boundary conditions uniquely determine the functions (j){r,v), g{r,v) and f{r,v). 

2.2 Structure of the equations of motion 



The nonzero equations of motion ( p.4|) consist of four nontrivial Einstein equations 



Err, Erv, Ejjy and Eii (where the index i runs over the d — 1 spatial directions) to- 
gether with the dilaton equation of motion. For the considerations that follow below, 
we will find it convenient to study the following linear combinations of equations 



El 


= g'^E^r 


El 


= g^'^E,, 


Eec 


= g^'^E^r 




d-l 


Ed 


= ^Eii 




i=l 


Ea, 





(2.10) 



Note that the equations El and El are constraint equations from the point of view 
of V evolution. 

It is possible to show that E^ and '^^'^^^"'^ both automatically vanish whenever 
El = El = Eij, = 0. This implies that this last set of three independent equations 
- supplemented by the condition that rE^c = at any one value of r - completely 
exhaust the dynamical content of (|2.4]) . As a consequence, in the rest of this paper 
we will only bother to solve the two constraint equations and the dilaton equation, 
but take care to simultaneously ensure that rE^c = at some value of r. It will often 
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prove useful to impose the last equation at arbitrarily large r. This choice makes 
the physical content of rEec = manifest; this is simply the equation of energy 
conservation in our system. 

2.2.1 Explicit form of the constraints and the dilaton equation 

With our choice of gauge and notation the dilaton equation takes the minimally 
coupled form 

dr {f-'gdr(P) + {f-'dr4>) + dr {f-'O,^) = (2.11) 

Appropriate linear combinations of the two constraint equations take the form 



(9,0)' 



/ (2.12) 



Note that the equations ( |2.12| ) (together with boundary conditions and the en- 
ergy conservation equation) permit the unique determination of f{r,Vo) and g{r,Vo) 
in terms of 0(r, fo) and (j){r,Vo) (where Vq is any particular time). It follows that / 
and g are not independent fields. A solution to the differential equation set ( p.llj ) 
and (p.l2|) is completely specified by the value of on a constant v slice (note that 



the equations are all first order in time derivatives, so on the slice is not part of 
the data of the problem) together with the boundary condition 0o(^^)- 

2.3 Explicit form of the energy conservation equation 

In this section we give an explicit form for the equation Egc = at large r. We 



specialize here to = 3 but see Appendix B.l for arbitrary d. Using the Graham 



Fefferman expansion to solve the equations of motion in a power series in - we find 



f{r,v)=r\ 1- ' .4 




o' 1 M ^4 



r2 V384 



'Jo) 




.(n.)^r^fi-2i^-^.0,l)) (2.13) 



^4 ' 



^^It turns out that both Ed and the dilaton equation of motion are automatically satisfied when- 
ever Eec together with the two Einstein constraint equations are satisfied. Consequently Eec plus 
the two Einstein constraint equations form another set of independent equations. This choice of 
equations has the advantage that it does not require the addition of any additional condition anal- 
ogous to energy conservation. However it turns out to be an inconvenient choice for implementing 
the e expansion of this paper, and we do not adopt it in this paper. 
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where the functions M{v) and L{v) are undetermined functions of time that are, how- 
ever, constrained by the energy conservation equation Eec, which takes the exphcit 
form 

M = 0o(^(0of -^-^00-) (2.14) 

In all the equations in this subsection and in the rest of the paper, the symbol P 
denotes the derivative of P with respect to our time coordinate v. Solving for M{v) 
we have 

Miv) = \jyt ((0o)' + \ (0o)' - 30oi:(t)) (2.15) 

16 



2.4 The metric and event horizon at leading order 



Later in this section we will solve the equations of motion (|2.11 ), ( 2.12| ) and ( p. 14 ) 



in an expansion in powers of e, the amplitude of the forcing function 0o(^)- In this 
subsection we simply state our result for the spacetime metric at leading order in e. 
We then proceed to compute the event horizon of our spacetime to leading order in 
e. We present the computation of the event horizon of our spacetime before actually 
justifying the computation of the spacetime itself for the following reason. In the 
subsections below we will aim to construct the spacetime that describes black hole 
formation only outside the event horizon. For this reason we will find it useful below 
to have a prior understanding of the location of the event horizon in the spacetimes 
that emerge out of perturbation theory. 



16 



We note parenthetically that (2.14) may be rewritten as 



To° = \^,^C (2.16) 
where the value C of the operator dual to the scalar field (f) and the stress tensor Tap are given by 
C = lim (dnct) + d^4>) 

r — >oo 

Ti:^ lim ( - (K - 2)6^ ^ 



r — >oo 



Where 

K!^ = Extrinsic curvature of the constant r surfaces, K = Kj^ 



(2.18) 

Qj^ — Einstein tensor evaluated on the induced metric of the constant r surfaces 



yielding 



T° = -2T^ = -2Ty = M{v) 

3 • (2-19) 

3L{v) - d^c^o 
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We will show below that to leading order in e, our spacetime metric takes the 
Vaidya form ( |1 . 1|) . The mass function M{v) that enters this Vaidya metric is also 
determined very simply. As we will show below, it turns out that L[v) ~ 0{e^) on 
our perturbative solution. It follows immediately from (|2.15 ) that the mass function 



M{v) that enters the Vaidya metric, is given to leading order by 

M{v) = C2{v) + 0{e^) 

1 r . ... (2.20) 

C2{v) = -- dtMt)<Po{t) 

(Here C2 is the approximation to the mass density, valid to second order in the 
amplitude expansion, see below). 

Note that, for v > 5t, C2{v) reduces to a constant M = C2 given by 

In the rest this subsection we proceed to compute the event horizon of the leading 
order spacetime ( |1.1| ) in an expansion in ea expansion. Let the event horizon manifold 
of our spacetime be given by the surface S = r — rniv) = 0. As the event horizon is 
a null manifold, it follows that dnSd^Sg^^ = 0, and we find 



dv 2 \ ^ff('^) 



(2.22) 



As M{v) reduces to the constant M = C2 ioi v > 6t, it follows that the event horizon 
must reduce to the surface rn = at late times. It is then easy to solve (|2.22|) for 
V < and v > 6t] we find 

rH{v) = M^, v>5t 

rniv) = Mh{^), < v < St 23) 
1 1 

-V^ ; , V<0 



ruKv) M3a;(0) 
where x{y) obeys the differential equation 



dx _ / M{y5t)\ 
dy~^Y\ Mx^ ) 



(2.24) 



a = M'i6t ~ es 



and must be solved subject to the final state conditions x = 1 for y = 1. (|2.24|) is 
easily solved in a perturbation series in a. We set 

x(y) = l + 5^«'^x„(y) (2.25) 
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and solve recursively for a;„(t). To second order we find^'^ 



1 /l _ M{z5t) 

^ M 



y 



2 



X2iy) = - I dz xi{z) ( 1 + 2^ 



y 



M{z6t)_ 



(2.26) 



In terms of which 

rniv) = (l + a xi{^^) + a^X2{-^) + 0{a^)^ {0 < v < St) (2.27) 

Note in particular that, to leading order, rniv) is simply given by the constant 
M5 for all v>0. 

2.5 Formal structure of the expansion in amplitudes 



In this subsection we will solve the equations ( ^.11| ), ( |2.12| ) and ( ^.141 ) in a pertur- 



bative expansion in the amplitude of the source function (j)o{v). In order to achieve 
this we formally replace 0o(^) with e0o(f) and solve all equations in a power series 
expansion in e. At the end of this procedure we can set the formal parameter e to 
unity. In other words e is a formal parameter that keeps track of the homogeneity of 
(pQ. Our perturbative expansion is really justified by the fact that the amplitude of 
00 is small. 

In order to proceed with our perturbative procedure, we set 



n=0 

oo 

g{r,v) = J2^''9n{r,v) (2.28) 

n=0 
oo 

(j){r,v) = ^e"0n(r,w) 



n=0 

with 

/o(r, v) = r, go{r, v) = 0o(r, v) = 0. (2.29) 

We then plug these expansions into the equations of motion, expand these equations 
in a power series in e, and proceed to solve these equations recursively, order by order 
in e. 



"'^^In this section we only construct the event horizon for the Vaidya metric. The actual metrics 
of interest to this paper receive corrections away from the Vaidya form, in powers of MSt. Conse- 
quently, the event horizons for the actual metrics determined in this paper will agree with those of 
this subsection only at leading order in MSt. The determination of the event horizon of the Vaidya 
metric at higher orders in MSt, is an academic exercise that we solve in this subsection largely 
because it illustrates the procedure one could adopt on the full metric. 



- 19 - 



The formal structure of this procedure is famihar. The coefficient of e" in the 
equations of motion take the schematic form 

H^^xUr,v)) = sl (2.30) 

Here xh stands for the three dimensional 'vector' of n*^ order unknowns, i.e. Xn — fn, 
Xn = 9n and Xn = The differential operator ifj is universal (in the sense that it 
is the same at all n) and has a simple interpretation; it is simply the operator that 
describes linearized fluctuations about AdS space. The source functions are linear 
combinations of products of x]n {fn < n) ; the sum over m over fields that appear in 
any particular term adds up to n. 

The equations ( p.30| ) are to be solved subject to the large r boundary conditions 



(2.31) 



lim 0i(r, v) = (f)o{r) 

Mr,v) <0{l/r), n>2 
fn{r,v)<Oil/r), n>l 
gn{r,v) < 0(r), n>l 

together with the initial conditions 

0n(r, v) = gn{r, v) = /,(r, v) = ioi v < {n > 1) (2.32) 

These boundary and initial conditions uniquely determine 0„, gn and /„ in terms of 
the source functions. 

All sources vanish at first order in perturbation theory (i.e the functions s\ are 
zero). Consequently, the functions /i and gi vanish but (pi is forced by its boundary 
condition to be nonzero. As we will see below, it is easy to explicitly solve for the 
function 0i. This solution, in turn, completely determines the source functions at 
C(e^) and so the equations ( ^.30] ) unambiguously determine g2, 02 and f2- This story 



repeats recursively. The solution to perturbation theory at order n — 1 determine 
the source functions at order n and so permits the determination of the unknown 
functions at order n. The final answer, at every order, is uniquely determined in 
terms of (f)o{v). 

To end this subsection, we note a simplifying aspect of our perturbation theory. 
It follows from the structure of the equations that 0„ is nonzero only when n is odd 
while fm and g^ are nonzero only when m is even. We will use this fact extensively 
below. 

2.6 Explicit results for naive perturbation theory to fifth order 

We have implemented the naive perturbative procedure described above to 0{e^). 
Before proceeding to a more structural discussion of the nature of the perturbative 
expansion, we pause here to record our explicit results. 
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At leading (first and second) order we find 



01 (r,^) = (l)o{v) + — 
r 



^ (2.33) 



C2{v) 3-2 



At the next order 



1 r 

M'f^^'t^) J B{x) dx 

^ 3^ ~ /I ^^""^ ^"^1 
,4(r,.) = ^ + ^{-0^ + 3 r B{x)dx 



^ 48r3 



^— {^B{v)^o - 40^00 + 300 ^ B{t)dt^ 



while 05 is given by 



<Pb{r,v) = I Bi{x) dx 



If" 5 f" n 



— 00 '/ — CxD '/ — 00 

^ /^^^(a;) dx + -^i I dy I B^ix) dx 

ry 



In the equations above 



(2.34) 



(2.35) 



B{v) 
Biiv) 



-C2{v) + 0000 



9 7 • •■ 

-^C2iv) + -0000 



B{x) dx 



+ lc2{v)4 + l^lB{v)-^4^o 
B2{v) = C4(t;)0o 

Bsiv) = ^ (-304 I B{x) dx + 74 



(2.36) 
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and the energy functions C2{v) and C^lv) (obtained by integrating the energy con- 
servation equation) are given by 

r 1 • ••• 

C^v) = J dt-<Po (-<Pl + j B{x) dx\ 

For use below, we note in particular that sA, v — 5t the mass of the black brane 
is given by C2{5t) — Ci{5t) + 0{e^) while the value of the dilaton field is given by 



1 r 

cl>{r,St)^—J B{x)dx 

+ d^J dyj B,{x)dx (2.38) 



' — OO d — oo 

5 n 1 z"'^* 

^24H / ^^(^)dx+— / Bz{x) dx 



1 /■'" 



2.7 The analytic structure of the naive perturbative expansion 

In this subsection we will explore the analytic structure of the naive perturbation ex- 
pansion in the variables v (for v > 5t) and r. It is possible to inductively demonstrate 
that 

• 1. The functions 02n+i, 5'2n+2 and /2n+2 have the following analytic structure 
in the variable r 

2n— 2 , iu / \ 



fc=0 ^ 



2n— 6 pj^ / \ 

hn{r,v)=rJ2¥^, in>S) (2.39) 



k=0 ^ 



(»>3) 

Moreover, when v > St (j)i{r, v) — /2(r, v) — fi{r, v) — Q while g2{T, v) — —^^^y^ 
and g4{r,v) = 

2. The functions 02n+i('^)' f2ni'^) 5'2n('^) functionals of (t)o{v) that 

scale like A^^^''^^'^, \-'^'^+^ and respectively under the scaling v — > A^;. 
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• 3. For V > 6t the functions 02n+i(^) polynomials in f of a degree that 

grows with n. In particular the degree of 02rt+i niost n — 1 + k; the degree 
of is at most n — 3 + k and the degree of gff^ is at most n — 4 + k. 

The reader may easily verify that all these properties hold for the explicit low 
order solutions of the previous subsection. 

2.8 Infrared divergences and their cure 

The fact that 02n+i('i^) arc polynomials in time whose degree grows with n immedi- 
ately implies that the naive perturbation theory of the previous subsection fails at 
late positive times. We pause to characterize this failure in more detail. As we have 
explained above, the field ^(r, v) schematically takes the form 



where (f>2n+i ~ ^(st)^^ at large times. Let us examine this sum in the vicinity r 
a surface that will turn out to be the event horizon of our solution. The term with 

2 2 

labels n, k scales like e x (e^ Now ^ = T is approximately the temperature 

of a black brane of event horizon rn- We conclude that the term with labels n, k 
scales like (vT)"~^+*^. It follows that, at least in the vicinity of the horizon, the 
naive expansion for is dominated by the smallest values of n and k when StT <^ 1. 
On the other hand, at times large compared to the inverse temperature, this sum is 
dominated by the largest values of k and n. As the sum over n runs to infinity, it 
follows that naive perturbation theory breaks down at time scales of order T~^. 

A long time or IR divergence in perturbation theory usually signals the fact 
that the perturbation expansion has been carried out about the wrong expansion 
point; i.e. the zero order 'guess' with which we started perturbation (empty AdS 
space) does not everywhere approximate the true solution even at arbitrarily small e. 
Recall that naive perturbation theory is perfectly satisfactory for times of order 5v 
so long as r » ^ . Consequently this perturbation theory may be used to check if our 
spacetime metric deviates significantly from the pure AdS in this range of r and at 
these early times. The answer is that it docs, even in the limit e — > 0. In order to see 
precisely how this comes about, note that the most singular term in g2n is of order 
r X for n > 1, the exact value of Qq = r"^ = {r x x r) . In other words Qq happens 
to be less singular, near r = 0, than one would expect from an extrapolation of the 
singularity structure of Qn at finite n down to n = 0. As a consequence, even though 
go is of lowest order in e, at small enough r it is dominated by the most singular 

2 

term in g2{r,v). Moreover this crossover in dominance occurs a-t ~ ^ ^ ^ and 
so occurs well within the domain of applicability of perturbation theory. In other 
words, in the variable range r ^ ^, g{r,v) is not uniformly well approximated by 




n,k 



2 
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5^0 = f'^ at small e but instead by 

(y((r, v) ^ r'^ — 

This implies that, in the appropriate parameter range, the true metric of the space- 



2 C,{V) 



time is everywhere well approximated by the Vaidya metric (|1 . 1| ) , with M{v) given 
by (12^201) in the limit e ^ 0. 

Of course even this corrected estimate for (yf(r, v) breaks down at r ~ ^. However, 
as we have indicated above, this will turn out to be irrelevant for our purposes as 

2 

our spacetime develops an event horizon at r ~ 

We will now proceed to argue that the metric is well approximated by the Vaidya 
form at all times (not just at early times) outside its event horizon, so that the Vaidya 
metric ( |1 . 1| ) rather than empty AdS space, constitutes the correct starting point for 
the perturbative expansion of our solution. 

2.9 The metric to leading order at all times 

The dilaton field and spacetime metric begin a new stage in their evolution at f = 
5t. At later times the solution is a normalizable, asymptotically AdS solution to 
the equations of motion. This late time motion is unforced and so is completely 
determined by two pieces of initial data; the mass density M[5t) and the dilaton 
function (j)[r,5t). As the naive perturbation expansion described in subsection 
is valid at times of order 5t, it determines both these quantities perturbatively in e. 
The explicit results for these quantities, to first two nontrivial orders in e, are listed 
in (PI). 

The leading order expression for the mass density is simply given by C2 in (|2.20| ) . 



Now if one could ignore 0(r, 5t) (i.e. if this function were zero) this initial condition 
would define a unique, simple subsequent solution to Einstein's equations; the uni- 
form black brane with mass density C2. While 0(r, 5t) is not zero, we will now show 
it induces only a small perturbation about the black brane background. 

In order to see this it is useful to move to a rescaled variable f = In terms 

of this rescaled variable, our solution aX v = 5t is & black brane of unit energy 
density, perturbed by 0(r, 5t). With this choice of variable the background metric is 
independent of e, so that all e dependence in our problem lies in the perturbation. 
It follows that, to leading order in e ( recall 0i(r, 5t) = 0) 

U) ^ ^ (1 + 0(.i)) 4x^(1^ Oi.i)) . ± (2.40) 



where, from subsection 2.6 



1 Z"^* 

4{6t) = - I B{x) dx (2.41) 



00 
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The important point here is that the perturbation is proportional to e and so 
represents a small deformation of the dilaton field about the unit energy density 
black brane initial condition. Moreover, any regular linearized perturbation about 
the black brane may be re expressed as a linear sum of quasinormal modes about 
the black brane and so decays exponentially over a time scale of order the inverse 
temperature. It follows that the initialy small dilaton perturbation remains small at 
all future times and in fact decays exponentially to zero over a finite time. The fact 
that perturbations about the Vaidya metric (|1 . 1|) are bounded both in amplitude as 
well as in temporal duration allows us to conclude that the event horizon of the true 
spacetime is well approximated by the event horizon of the Vaidya metric at small 
e, as described in subsection 0|. 



2.10 Resummed versus naive perturbation theory 

Let us define a resummed perturbation theory which uses the corrected metric ( |l.lj ) 
(rather than the unperturbed AdS metric) as the starting point of an amplitude 
expansion. This amounts to correcting the naive perturbative expansion by working 
to all orders in M ~ e^, while working perturbatively in all other sources of e depen- 
dence. As we have argued above, resummed perturbation theory (unlike its naive 
counterpart) is valid at all times. 

We have seen above that the naive perturbation theory gives reliable results when 
vT <C 1. This fact has a simple 'explanation'; we will now argue that the resummed 
perturbation theory (which is always reliable at small e) agrees qualitatively with 
naive perturbation theory vT <^ 1. 

At each order, resummed perturbation theory involves solving the equation 



dr 



+ 2rdvdr{r(f)) = source (2.42) 



The naive perturbation procedure requires us to solve an equation of the same form 
but with M set to zero. In the vicinity of the horizon, the two terms in the expression 
(1 — ^^:^) are comparable, so that the resummed and naive perturbative expansions 
can agree only when the entire first term on the LHS of ( |2.42| ) is negligible compared 
to the second term on the LHS of the same equation. The ratio of the first term 
to the second may be approximated by rv where v is the time scale for the process 
in question. Now the term multiplying the mass in ( p.42| ) is only important in the 



neighborhood of the horizon, where r ~ Ma ~ T where T is the temperature of the 



^^This is conceptually similar to the coupling constant expansion in finite temperature weak 
coupling QED. There, as in our situation, naive perturbation theory leads to IR divergences, which 
are cured upon exactly accounting for the photon mass (which is of order gyAi)- Resummed 
perturbation theory in that context corresponds to working with a modified propagator which 
effectively includes all order effects in the photon mass, while working perturbatively in all other 
sources of the fine structure constant a. 
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black brane. It follows that resummed and naive perturbation expansions will differ 
substantially from each other only at time scales of order and larger than the inverse 
temperature. 

Let us restate the point in a less technical manner. The evolution of a field 0, 
outside the horizon of a black brane of temperature T, is not very different from 
the evolution of the same field in Poincare patch AdS space, over time scales v 
where vT <^ 1. However the two motions differ significantly over time scales of 
order or greater than the inverse temperature. In particular, in the background of 
the black brane, the field outside the horizon decays exponentially with time over 
a time scale set by the inverse temperature; i.e. the solution involves factors like 
e~^^. As the temperature is itself of order es, naive perturbation theory deals with 
these exponentials by power expanding them. Truncating to any finite order then 
gives apparently divergent behavior at large times. Resummed perturbation theory 
makes it apparent that these divergences actually resum into completely convergent, 
decaying, exponentials. 

2.11 Resummed perturbation theory at third order 

In the previous subsection we have presented explicit results for the behavior of the 
dilaton and metric fields, at small e and for early times vM^ <^ 1. The resummed 
perturbation theory outlined in this section may be used to systematically correct 

____ 2 

the leading order spacetime ( |1.1D at all times, in a power series in es. In this section 
we explicitly evaluate the leading order correction in terms of a universal (i.e. 0o 
independent) function ip^Xjy), whose explicit form we are able to determine only 
numerically. 

Let us define the function ipi^j v) as the unique solution of the differential equa- 
tion 

(x^ (l - ^) d,i) + 2xdyd, (xij) = (2.43) 



subject to the boundary condition ~ ^{^) large x and the initial condition 
^jJ{x,0) = ^. The leading order solution to the resummed perturbation theory for 
0, for t> > is given by 

ivi M 3 



Unfortunately, the linear differential equation ( p. 43 ) - appears to be difficult to 



solve analytically. In this section we present a numerical solution of ( p.43|) . Although 



we are forced to resort to numerics to determine %p{x,y), we emphasize that a single 
numerical evaluation suffices to determine the leading order solution at all values 
of the forcing function (I)q{v). This may be contrasted with an ab initio numerical 
approach to the full nonlinear differential equations, which require the re running 
of the full numerical code for every initial function 0o- In particular the ab initio 
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numerical method cannot be used to prove general statements about a wide class of 
forcing functions 0o- 

In Figure || we present a plot of ip{^,y) against the variables u and y. The 
exterior of the event horizon lives in the compact interval ^ = m G (0, 1), and in our 
figure y runs from zero to three. 




Figure 3: Numerical solution for dilaton to the leading order in amplitude at late time 



In order to obtain this plot we rewrote the differential equation (|2.43| ) in terms of 
the variable u = ^ (as explained above) and worked with the field variable y) = 
{1 — u)iIj{^, y). Recall that our original field ip is expected to be regular at the horizon 
M = 1 at all times. This expectation imposes the boundary condition x(. 999999, y) = 
0. We further imposed the condition of normalizability x(0,?/) = and the initial 
condition x{u, 0) = (0.999999 — u)u^. Of course 0.999999 above is simply a good 
approximation to 1 that avoids numerical difficulties at unity. The partial differential 
equation solving routine of Mathematica-6 was able to solve our equation subject to 
these boundary and initial conditions, with a step size of 0.0005 and an accuracy goal 
of 0.001; we have displayed this Mathematica output in figure |[ In order to give a 
better feeling for the function ip{x, y) in figure^ we present a graph of '(/'((fy, y) (i-e. as 
a function of time at a fixed radial location). Notice that this graph decays, roughly 
exponentially for v > 0.5 and that this exponential decay is dressed with a sinusodial 
osciallation, as expected for quasinormal type behavior. A very very rough estimate 



-27- 



^/r(1.42857, y) 




y 



Figure 4: A plot of ilj{jjj,y) as a function of y 



of this decay constant uj may be obtained from equation 



^(0^7,-5) 



e which gives 



uj ^ 8.9T (here T is the temperature of our black brane given by T = ^). This 
number is the same ballpark as the decay constant for the first quasi normal mode 
of the uniform black brane, uj = 11.16T, quoted in [H0|. 



3. Spherically symmetric asymptotically flat collapse 
3.1 The Set Up 

In this section^^ we study spherically symmetric asymptotically flat solutions to Ein- 
stein gravity (with no cosmological constant) interacting with a minimally coupled 
massless scalar field, in 4 bulk dimensions. The Lagrangian for our system is 



S = I d'x^(^R-^{d<pY^ (3.1) 



We choose a gauge so that our metric and dilaton take the form 



ds = 2drdv — g{r,v)dv +f {r,v)dQ, 



2 



(3.2) 



where dfl2 is the line element on a unit two sphere. We will explore solutions to the 



0(r,w). 
t on a i; 

equations of motion of this system subject to the pure flat space initial conditions 

(3.3) 



g{r,v) 


= 1, 


{v < 0) 


f{r,v) 


= r, 


{v < 0) 


(j){r,v) 


= 0, 


{v < 0) 



-^^We thank B. Kol and O. Aharony for discussions that led us to separately study collapse in fiat 
space. 
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and the large r boundary conditions 



g{r,v) = l + 0{-) 
r 

S{r,v) = r(l + 0(^)\ (3.4) 
^(r..) = ^ + 0(J,) 

where il){v) takes the form 

^{y) = 0, [V < 0) 

i){v) < ef5t, {0<v < 5t) (3.5) 

^(t;) = {V> 5t), 

In other words our spacetime starts out in its vacuum, but has a massless pulse of 
limited duration focused to converge at the origin at f = 0. This pulse could lead to 
interesting behavior - like black hole formation, as we explore in this section. 

The structure of the equations of motion of our system was described in subsec- 
tion p^ . As in that subsection, the independent dynamical equations for our system 
may be chosen to be the dilaton equation of motion plus the two constraint equa- 
tions, supplemented by an energy conservation equation. The explicit form of the 
dilaton and constraint equations is given by 

Or [fgdrCf) + {fdA) + dr {fd,<P) = 

{dAf = (3.6) 

dr ifgdrf + 2fdJ) = 1 

As in the previous section, we may choose to evaluate the energy conservation equa- 
tion at large r. As we have explained, the large r behavior of the function g is given 

by 

gir,v) = l-^ + 0{^^) (3.7) 
The energy conservation equation, evaluated at large r, yields 

M = -^ (3.8) 

The equations ( pl6|) together with ( |3.8|) constitute the full set of dynamical equations 
for our problem. 

By integrating (|3.8| ) we find an exact expression for M{v) 

M{v) = (3.9) 
Note in particular that M[v) reduces to a constant M for v > 5t where 

M = ^-^^ ~ e}6t (3.10) 
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3.2 Regular Amplitude Expansion 

Our equations may be solved in the amplitude expansion formally described in 
i.e. in an expansion in powers of the function ip{v). As we will argue in this paper, 
there are two inequivalent valid amplitude expansions of these equations. In the 
first, the spacetime is everywhere regular and the dilaton is everywhere small. In 
the second, the spacetime is singular at small r but this singularity is shielded from 
asymptotic infinity by a regular event horizon. The second amplitude expansion 
reliably describes the spacetime only outside the event horizon; this expansion works 
because the dilaton is uniformly small outside the event horizon. As we will see two 
amplitude expansions described above have non overlapping regimes of validity, and 
so describe dynamics in different regimes of parameter space. 

In this subsection we briefly comment on the more straightforward fully regular 
expansion. At every order in perturbation theory, the requirement or regularity 
uniquely determines the solution. Explicitly at first order we have 

Mr,v) = (3.11) 

r 

The perturbation expansion that starts with this solution is valid only when 0(r) is 
everywhere small. 0(r) reaches its maximum value near the origin, and 0i(O,f) ~ 
2ip{v) ~ ef. Consequently the regular perturbation expansion, sketched in this 
section, is valid only when <^ 1 i.e. when ^ ^ 1. 
At next order in the amplitude expansion we find 



1 



f2{r,v) = -yrj p[dp(j)i{p,v)] dp - J P [dpMp^'^)] dp 
92{r,v) = -2dj2{r,v) 5,./2(r,f) 



(3.12) 



The integration limits in the expression for f2{r, v) in p.l2| are fixed such that at large 
r /(r, v) decays like K The integration constant in (72 (r, v) is fixed by the requirement 
that the solution be regular at r = 0. 

3.2.1 Regularity implies energy conservation 

In this subsection we pause to explain an interesting technical subtlety that arises in 
carrying out the regular amplitude expansion. The discussion of this subsection will 
play no role in the analysis of spacetimes that describe black hole formation, so the 
reader who happens to be uninterested in the regular expansion could skip to the 
next section. 

Note that in order to obtain ( |3.12| ) we did not make any use of the energy 
conservation equation. We will now verify (first in terms of the answer, and then 
more abstractly) that ( |3.12| ) automatically obeys the energy conservation equation. 
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At large r, these functions have the following expansion 



0i(r,f) 



8r 
C2{v) 



(3.13) 



where 



f2i0,v) 



If our solution does indeed obey the energy conservation relation, we should find that 
C2{v) is equal to M{v) in ( 3.10| ). We will now proceed to directly verify that this is 
the case. 

The first term in C2(f) comes from the coefficient of ^ in (9„/2(r, f). For the 
second term in the expression for C2{v), /2(0), is given by 



1 /"°° 

/2(0,w) = -4 y [dpMp,'v)f dp 



The integrand in this expression may be split into four terms in the following way. 



[dr(f)i{r,v)f = 2il){v)dr 



2ip{v)dr 



iIj{v — 2r) 



iIj{v — 2r) 



+ 



a, 



il){v — 2r) 



1 2 



+ 4 



iIj{v — 2r) 



— dr 



ip'^{v — 2r) 



(3.14) 



Now each of the terms can be integrated. 



\p{v~2p) 



P 



dp 



-21im 



'il){v)%l){v — 2r) 



-21im 



r^O T 



P 

4 



dp = lim 







2 



^Ij{v - 2p)J dp = 2 I i>{tf dt 
^^{v-2p) 



(3.15) 



P 



i'(v-2rY ?/'(f)^ 
dp = lim = lim 



Adding all the terms one finally finds 



-f2{0,v) 



1 



This implies 



^/'(t) dt 



?/'(t) dt = M(f ) 



(3.16) 



(3.17) 
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as expected from energy conservation. 

Let us summarize In order to obtain our result for g2 above, we were required to 
fix the value of an integration constant. The value of this constant may determined 
in two equally valid ways 

• By imposing the energy conservation equation Ef.c 

• By demanding regularity of the solution at r = 

In fact these two conditions are secretly the same, as we now argue. As we have 
explained in subsection p.2| , dr{rEec) automatically vanishes whenever the three 
equations ( p.6|) are obeyed. Consequently, if rEec vanishes at any one value of r 
it automatically vanishes at every r. Now the equation E^c evaluates to a finite value 
at r = provided our solution is regular at r = 0. It follows that the regular solution 
automatically has rEec = everywhere. 

Configurations in the amplitude expansion of the previous section (or the sin- 
gular amplitude expansion we will describe shortly below), on the other hand, are 
all singular at r = 0. rE^c does not automatically vanish on these solutions, and 
the energy conservation equation E^c is not automatic but must be imposed as an 
additional constraint on solutions. 

It would be a straightforward - if cumbersome - exercise to explicitly implement 
the perturbation theory, described in this subsection, to higher orders in e/. As our 
main interest is black hole formation, we do not pause to do that. 

3.3 Leading order metric and event horizon for black hole formation 

In the rest of this section we will describe the formation of black holes in flat space 
in an amphtude expansion. In contrast with the previous subsection, our amplitude 
expansion will be justified by the small parameter ^. Our analysis will reveal that 
our spacetime takes the Vaidya form to leading order in 4-, 

ds^ = 2drdv - (^1 - ^^^^ dv^ + r^dQl (3.18) 

where M{v) is given by (|3.9|) . 

In this subsection we will compute the event horizon of the spacetime ( p.l^ ) 
at large e/. We present the computation of this event horizon even before we have 
justified the form (|3.18|) , as our aim in subsequent subsections is to have a good per- 
turbative expansion of the true solution only outside the event horizon; consequently 
the results of this subsection will guide the construction of the amphtude expansion 
in subsequent subsections. 



- 32 - 



As in the previous section the event horizon takes the form 
rniv) = M, {v> 6t) 

rniv) = Mx{^), {0 < v < 6t) (3.19) 
ot 

rniv) = Mx{0) + v, (-a;(0) < v < 0) 

where the function x(t) may easily be evaluated in a power series in |j ~ 4-- We 
find 

(3.20) 




In particular th = M for all > at leading order. 

3.4 Amplitude expansion for black hole formation 

Let us now construct an amplitude expansion (i.e. expansion in powers of ip{v)) 
of our solution in the opposite limit to that of the previous subsection, namely 
~ ej ^ 1. It is intuitively clear that such a dilaton shell will propagate into 
its own Schwarzschild radius and then cannot expand back out to infinity. In other 
words the second term in (|3.11|) cannot form a good approximation to the leading 



order solution for the collapse of such a shell. Now (|3.11|) deviates from 



0i(r,t;) = ^; (3.21) 



only at spacetime points that feel the back scattered expanding wave in (|3.11| ). This 



observation suggests that ( p.21|) itself is the appropriate starting point for the am- 



plitude expansion at large e/, and this is indeed the case. 



The incident dilaton pulse ( |3.21| ) will back react on the metric; above we have 
derived an exact expression for one term - roughly the Newtonian potential - (see 
( |3.10| )) of this back reacted metric. Including this backreaction (all others turn out 
to be negligible at large e/) the spacetime metric takes the form 

ds^ = 2dvdr - dv\l - ^^) + r^dnl (3.22) 



As we have explained in the previous subsection, this solution has an event horizon 
located at r^^ ~ M ~ e'jSt for i; > (see below). Consequently, 0i(r, f) outside 

the event horizon < ~ — ~ \/—, i.e. is parametrically small at large et. This 
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fact allows us to construct a large ef amplitude expansion for the solution outside 
its event horizon. 

The perturbation expansion of our solutions in || is similar in many ways to 
the perturbation theory described in detail in section 0. As in that section, the true 
(resummed) expansion (built around the starting metric ( p.22| )) is well approximated 
at early times by a naive expansion built around unperturbed flat space. Naive and 
resummed expansions agree whenever the first term in the first equation of ( |3.(j| ) 
is negligible compared to the other terms in that equation, i.e. for f <C M ~ 
^"fSt. As e/ is large in this subsection, naive and resummed perturbation theory are 
simultaneously valid for times that are of order 6t. However we expect the naive 
expansion to break down at v ^ M. We will now study the naive expansion in more 
detail and confirm these expectations. 

3.5 Analytic structure of the naive perturbation expansion 

In this subsection we describe the structure of a perturbative expansion built starting 
from the flat space metric. We expand the full solution as 

oo 

0(r, v) = J2 "^Sn+l 

n=0 

oo 

f{r,v) = r + J2F2n{r,v) (3.23) 

n=l 

oo 

g{r,v) = 1 + ^G2n{r,v) 

n=l 

where, by definition, the functions $rra Fm and Gm are each of homogeneity m in the 
source function ip{v). As explained above we take 

$i(r,t;) = ^ (3.24) 
r 

By studying the formal structure of the perturbation expansion, it is not difficult to 
inductively estabhsh that 

• 1. The functions $2n+i, and G2n have the following analytic structure in 
the variable r 



m=0 
oo 



F2n{r.-) = rj:^ (3.25) 



r 

m=0 



771=0 
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• 2. The functions ^2n+ii 

(f ) are each functionals of ip{v) that 
scale hke A"* A™ and A"*"^ under the the scahng v —>■ Xv. 



• 3. For V > 6t the $2n+i(^) polynomials in v of degree < n + m — 1; F2n{v) 
and are polynomials in v of degree < n + m — 3 and n + m — 4 respectively. 



It follows that, say, 4>{r,v), is given by a double sum 



n n,m=0 



Now sums over m and n are controlled by the effective expansion parameters ~ ^ 
(for m) and ~ ~ (for n; recall that in the neighborhood of the horizon 
th ~ Stej). 

It follows that the sum over ni is well approximated by its first few terms if only 
V <^ M (recall we are interested in the solution only for r > M). The sum over 
n may also be truncated to leading order only for v <^ M. As anticipated above, 
therefore, our naive perturbation expansion breaks over time scales v of order and 
larger than M. 

Let us now focus on times v of order St. Over these time scales naive perturbation 
theory is valid for r ^ efSt (recall that this domain of validity includes the event 
horizon surface rn ~ elSt). Focusing on the region of interest, r ^ r^, sn+m+i 

scales 

like 2n+2m+i ■ It follows that . i, with equal values of n + m are comparable at 
times of order 5t. For this reason we find it useful to define the resummed fields 

02n+l(r, V)=}_^ _^ 



k=0 ^ 



n-2 pfe 

I 

k=0 



hn{r,v)=:rS„,FS + r'Y:^^§^ (3.26) 

fc=0 

n-2 ^fc 



<p2n-ii unlike $2n-i; TGceives contributions from only a finite number of terms at any 
fixed n, and so is effectively computable at low orders. According to our definitions, 
0m, fm and Qm Capture all contributions to our solutions of order at time scales 
of order 5t. 

We now present explicit computations of the fields 0^, fm and Qm up to 5th 
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order. We find 



Mr,v) 



8r 
M{v) 

r 



384r^ 32r^ 
, , ip(v)ilj(vf M(v)ih(vV ip(v)B(v) 

^ ( ^ ^^^^ 



(48i?(t) - 16^(t)3) rft 



192r 



4 



^oo ^W^W + 21i?W} + 3M(t) {^(t)3 - 18i?(t)} 



192r5 



= f ^{t) (^-M{t)+^{t)^{t)^ dt 



+ 

Where 



3.6 Resummed perturbation theory at third order 

As in the previous subsection, even at times of order 6t (where naive perturbation 
theory is valid) naive perturbation theory yields a spacetime metric that is not uni- 
formly well approximated by empty flat space over its region of validity r ^ ^t^f- 
The technical reason for this fact is very similar to that outlined in the previous 
section; go is a constant, so is smaller at r ~ r^/ than one would have guessed from 
the naive extrapolation of ( |3.25| ) to n = 0. It follows that, in the previous section 
that, even at arbitrarily small e, the resultant solution is well approximated by 

g{r, V) ^ 1 — 



rather than the flat space result g{r, v) = 1, over the full domain of the amplitude 
expansion. It follows that the correct (resummed) amplitude expansion should start 
with the Vaidya solution ( |3.22| ) rather than the empty flat space. The IR divergences 



of the naive expansion are a consequence of the incorrect choice of starting point for 
the perturbative expansion. 

At V = 6t our metric, to leading order, is the Schwarzschild metric of a black 
hole Schwarzschild radius M with a superposed dilaton (and consequently metric) 
perturbation. We will now demonstrate that these pertubrations are small. As in 
the previous section, it is useful to define rescaled radial and time variables x = 
and y = j^- In terms of the rescaled variables, the leading order metric takes the 
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form 



ds^ = (^dxdy - dy^ (^1 - + x^dQj^ (3.28 
while the (p perturbation is given to leading order by 

<pm _ <pi{6t) 1 



(3.29) 



(recall from ( 3.27|) that 



0O(5t) = ] I ij{v) -M{v) + ij{v)ij{v) dv (3.30) 



-St 

and M(f) is given in ( |3.9| )). 

As a constant rescaling of the metric is an invariance of the equations of motion 
of the Einstein dilaton system, the factor of in ( |3.28| ) is irrelevant for dynamics. 
As the dilaton perturbation above is parametrically small (0(l/ej)) the subsequent 
evolution of the dilaton field is linear to leading order in the ^ expansion. 

Let x(a^; y) denote the unique solution to 

(^x' (^1 - d,x^ + 2xdyd, (xx) = (3.31) 

subject to the boundary condition x ~ ^(^a) large x and the initial condition 
x(a;, 0) = The leading order solution to the resummed perturbation theory for 0, 
for V > 6t, is given by 

.-mjril^) (3.32) 



M3 ^ V^' M 

Unfortunately, the function x(x, y) appears to be difficult to determine analyti- 
cally. As in section |^ this solution may presumably be determined numerically with 
a little effort. We will not attempt the requisite numerical calculation here. In the 
rest of this subsection we will explain in an example how the general analysis of this 
subsection yields useful precise information about the subleading solution even in 
the absence of detailed knowledge of the function x(a;,?/). 

Consider a spherically symmetric shell, of the form discussed in this section, 
imploding inwards to form a black hole. On general grounds we expect some of 
the energy of the incident shell to make up the mass of the black hole, while the 
remaining energy is reflected back out in the form of an outgoing wave that reaches 

Let the fraction of the mass that is reflected out to X+ be denoted by /. / is 
one of the most interesting and easily measured observables that characterize black 
hole formation. 



^°In fancy parlance / = where A is the ADM mass of the spacetime and B is the late time 
Bondi mass. 
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At leading order in the expansion in ^ our spacetime metric takes the Vaidya 
form with no outgoing wave, and so / = 0. This prediction is corrected at first 
subleading order, as we now explain. It follows on general grounds that, at late 
times 

X{x,y) ^ 

X 

for some function C('^)- Note that (, like the function x, is universal (i.e. independent 
of the initial condition iplv)). It follows that at late times (and to leading order) 

= M^^iaa. (3.33) 



It then follows from (|3.1CI|) (but now applied to an outgoing rather than an ingoing 



wave) that the energy^^ carried by this pulse is 



It follows that 

f = A 



m 

M3 



oo 

2 



(3.35) 



oo 



( |3.35| ) analytically determines the dependence of / on the shape of the incident wave 
packet, tp{v) (recall that (p^iSt) and M are determined in terms of ip{v) by (^ 



and (|3.10|) ). Detailed knowledge of function x(a;, y) is required only to determine the 



precise value of universal dimensionless number A. 



4. Spherically symmetric collapse in global AdS 

We now turn to the study of black hole formation induced by an ingoing spherically 
symmetric dilaton pulse in an asymptotically AdSd+i space in global coordinates. 
As in section ^ our bulk dynamics is described by the Einstein Lagrangian with 
a negative cosmological constant and a minimally coupled dilaton. However as in 
section ^we study solutions that preserve an SO{d) invariance; this SO{d) may be 
thought of as the group of rotations of the boundary S'^~^. As in both sections ^ and 
^ our symmetry requirement determines our metric up to three unknown functions 
of the two variables; the time coordinate v and the radial coordinate r. Our solutions 
are completely determined by the boundary value, (f)o{v) of the dilaton field. As in 

^^We have chosen our units of energy so that a black hole with horizon radius rj^ has energy M . 
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section || we assume that (j)Q{v) is everywhere bounded by e and vanishes outside the 
interval (0, 6t). Through out this section we will focus on the regime St <^ R (where 
R is the radius of the boundary sphere) and e -C 1. The complementary regime 
5t ^ R and arbitrary e is under independent current investigation [[79[] . 

The collapse process studied in this section depends crucially on two independent 
dynamical parameters; x = ^ together with e of previous subsections. We study the 
evolution of our systems in a limit in which x and e are both small. The problem 
of asymptotically AdS spherically symmetric collapse is dynamically richer than the 
collapse scenarios studied in sections ^ and ^ and indeed reduces to those two special 
cases in appropriate limits. 

4.1 Set up and equations 

The equations of motion for our system are given by ( ^^.41 ). The form of our metric 
and dilaton is a slight modification of (|2.6|) 



ds^ = 2drdv - g{r, v)dv^ + f\r, v)dnl_, 
(p = 0(r, v). 

where dQ'^_-^ represents the metric of a unit d — 1 sphere. Our fields are subject to 
the pure global AdS initial conditions 

g[r,v) = r^ + {v<0) 

f{r,v)=rR, {v<0) (4-2) 
(j){r,v) = 0, {v<0) 

and the large r boundary conditions 

g{r,v) = r' (^1 + 
fir,v) = r(^R + Oi^)^ (4-3) 

r 



Equations ( |2.4| ), ( [4.1|) , ( |4.2| ) and ( |4.3D together constitute a completely well defined 
dynamical system. Given a particular forcing function (poi^), these equations and 
boundary conditions uniquely determine the functions (/>(r, f), g{r,v) and f{r,v). 

The structure of the equations of motion of our system was described in sub- 
section p?2| . In particular, we may choose the dilaton equation of motion, together 
with the two constraint equations, as our independent equations of motion; this set 
is supplemented by the energy conservation relation. With our choice of gauge and 
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notation, the dilaton equation of motion and constraint equations take the exphcit 
form 

dr [f-^gdA) + dv {f-'dr<P) + dr {f-'gd,<p) = 

(^.0f + ^i^^^O (4.4) 

dr {f-'gdrf + 2f-'dj) - d f-' -{d- 2)f-' = 

As in section 2, the initial data needed to specify a solution to these equations is 
given by the value of 0(r) on a given time slice, supplemented by the initial value of 
the mass, and boundary conditions at infinity. In order to obtain an explicit form 
for the energy conservation equation we specialize to c? = 3 and explicitly 'solve' our 
system at large r a la Graham and Fefferman. We find 



/(r,.;)=i?r|l-^ + 0(l; 

3^1 



j_ l_^_MW^„,i> (4.5) 



9(r,v) = r'{— + + 0( 



4' 



The energy conservation equation constrains the (otherwise arbitrary) functions 
M{v) and L{v) to obey 

M = (^UL{v) +4^-3 (0o)' + 40o(^) j (4-6) 

22 

4.2 Regular small amplitude expansion 

As in section ^ there are two legitimate amplitude expansions of spacetime we wish 
to determine. In this subsection we discuss the expansion analogous to the expansion 
of subsection p.2| . That is we expand all our fields as in ( |2.28| ) (where the functions 



^^Notc that the stress tensor and Lagrangian C of our system are given by 

t: = Miv) 

Te-T^- — (4.7) 
3 / • \3 



C = -3L{v)~^ + -[^o) ~'4>,{v) 



It follows that ( |4.6D may be rewritten as M 



<t>oC- 
2 ■ 
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fn, Qn and (pn are all defined to be of homogeneity n in the boundary field ^o) and 
demand that all functions are everywhere regular. The requirement of regularity, 
together with our boundary and initial conditions, uniquely specifies all functions in 
( p.28|) . Explicitly, to second order we find 



where 



m=0 



R 



^o{v - mnR) + — + 0o(^^ - Rmn - 2Rtan-^ (rR)) 



boiv - mnR - 2Rtan-^ (rR)) 



r / pK{p,v)dp 



1 

4r 



2r 



p^K{p, v) dp 



fOO 

pK{p,v) dp + 2r^ I p^K{p,v)dp 



+ 1^ p\^ + p')K{p,v)dp 



K{p,v) = {drMr,v)y 



The perturbation expansion in this section is valid only if 0(r, v) is everywhere 
small on the solution. (j)i{r,v) reaches its maximum value in the neighborhood of 
the origin where it is given approximately by 0o + 0o ~ e + Consequently the 
validity of the amplitude expansion sketched in this section requires both that e ^ 1 
and that x"^ ^ e. 

We have chosen integration constants to ensure that the solution in ( [4.8|) is 
regular at r = 0. In particular 



92(0, v) 



As in subsection ^.2| , this choice automatically implies the energy conservation equa- 
tion. In particular, expanding g2{r,v) at large r we find 



-M{v) = - 



4 



OO -| 

2/ I „2 



p(^ + p)^(p^ ^) - Mvy - 2Mv)Mv) 



(4.9) 



(this equation is valid only for v < tcR; it turns out that M{v) is constant for v > 6t) 
in agreement with the energy conservation equation. 

Finally, let us focus on the coordinate range rR, <^ 1 and also require that 
X is small so that the time scale in 0o is also smaller than the AdS radius. In this 
parameter and coordinate range ( [4.8| ) should reduce to a solution of the fiat space 
propagation equation ( ^.11| ); this is easily verified to be the case. In the given variable 
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and parameter regime, all terms with (^4.8|) with m 7^ vanish; tan ^{Rr) ^ rR and 
the first and the third terms in ( [4.8|) are negligible compared to the second and 



fourth as x is small. Putting all this together, (|4.8| ) reduces to (|3.11| ) under the 
identification ipi"^) = R'^4'o{v), once we also identify the coordinate r of section ^ 
with R^r in this section. Notice that this replacement implies that e/ = ^ (where 
tf was the perturbative expansion of section |^). This identification of parameters 
is consistent with the fact that the expansion of this section breaks down when -% 



becomes large, while the expansion of subsection |3.2| breaks down at large e/ 



4.3 Spacetime and event horizon for black hole formation 

In the rest of this section we will describe the process of black hole formation via col- 
lapse in an amplitude expansion. As in earlier sections, the spacetime that describes 
this collapse process will turn out to be given, to leading order, by the Vaidya form 



ds' = 2drdv -{^ + r^- 1 dv' + R'r'd^i, 

R^ r 

{r,v) = (f)o{v) + — 
r 



(4.10) 



where M(t>) is approximated by C2{v), the order piece of ( [4.6|) 

C2{v) = -\ £^ dtMt) (^0o(^) + ^ j (4.11) 

In this subsection we will compute the event horizon of the spacetime ( |4.10|) in a 
perturbation expansion in a small parameter, whose nature we describe below. The 
horizon is determined by the differential equation 

1^ MM 

dv ii^ Vff 

where M{v) reduces to a constant M for t > 6t. At late times the event horizon 
surface must reduce to the largest real solution of the equation 

1 . n.o M 



H 



It then follows from ( |4.12|) that 



^niv) = r%, (f > 5t) 

1) 

rniv) = r%x{—), (0 < 
tan"^ (rHiv)) = tan^^ {r%x{0)) +v {v < 0), tan-^(rH(f )) > 



rniv) = r%x{^), {0 < v < 6t) (4.13) 
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As in previous subsections, the function x{t) is easily generated in a perturbation 
expansion 

xit) = 1 + i^-^j x,{t) + (^^j x,{t) + ... (4.14) 
The small parameter for this expansion is . This parameter varies from approx- 

2 2 4 2 

imately es when x ^ es to ^ when x ^ es and is always small provided x -C ^/e 
and e <^ 1. These conditions will always be met in our amplitude constructions 
below. Note that the event horizon of our solution is created (at r = 0) at the time 
V = — tan~^(r^)+ subleading. 

Explicitly working out the perturbation series we find 



,1 _ MjySt) 
X^{t) = - I dt 



(4.15) 



4.4 Amplitude expansion for black hole formation 

The amplitude expansion of the previous subsection breaks down for x"^ <C e. As in 
section ^, we have a new amplitude expansion in this regime. As in section ^, the 
starting point for this expansion is the Vaidya metric and dilaton field (|4.10|) . As in 
sections || and |, the perturbation expansion based on (|4.10|) is technically difficult 



to implement at late times. However as in earlier sections, at early times - i.e. times 
of order 6t - the perturbative expansion is well approximated by the naive expansion 
based on the solution ( |4.10|) with M{v) set equal to zero. Following the terminology 
of previous sections we refer to this simplified expansion as the naive expansion. In 
the rest of this subsection we will elaborate on the analytic structure of the naive 
perturbative expansion. 

In order to build the naive expansion, we expand the fields f{r,v), g{r,v) and 
4>{r,v) in the form ( p.28| ). It is not too difficult to inductively demonstrate that 

• 1. The functions 02n+i, g2n and f2n have the following analytic structure in the 
variable r 

oo ^ 2n+m-2 ik,m i \ 

02n+l(r,t;) = 2^-^ 2^ (^>1) 
m=0 fc=0 

g2n[r,V) - +r 2^ ^2n-k+m l"- ^ 

m=0 k=0 
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• 2. The functions 4>2n+i{'^)^ /2n™('^) g^^iy) are functionals of (^^[v) that 
scale hke _\-2n.-i+m+fe^ y^-2n+m,+k j)^-2n+m+fc-i respectively under the scaling 

V —>■ Xv. 

• 3. For V > 5t we have some additional simplifications in structure. At these 
times /4(r, v) = and (74 (r, f ) = Further, the sums over k in the second 
and third of the equations above run from to 2n — 6 + m and 2n — 5 + m 
respectively. Finally, functions 02n+i('^) ^ire all polynomials in t> of a degree 
that grows with n. In particular the degree of 02n+i most n — 1 + k + m; 
the degree of is at most n — 3 + k + m and the degree of is at most 
n — 4 + k + m. 

As we have explained above, 

00 CO , 2n— 2+m ,fe,m / \ 
V-V- 1 V- 02;+lW 



n=l m=0 k=0 



We will now discuss the relative orders of magnitude of different terms in this summa- 
tion. Abstractly, at times that are larger than or of order St, the effective weighting 

2 

factor for the sum over n,m,k respectively are approximately given by ■:pij^) -j^ 
and vr respectively. We will try to understand the implications of these estimates in 
more detail. 

2 

Let us first suppose that x <^ es. In this case the black hole that is formed has 

2 

a horizon radius of order ^ ^ ^ { this estimate is corrected in a power series in 
%). Consequently, the resultant black hole is large compared to the AdS radius. At 
m = 0, this regime, the summation over k and n simply reproduce the solution of 
section ^. As in section ^ these summations are dominated by the smallest values 
of k and n for vr^ ~ vT <C 1, in the neighborhood of the horizon. As in section 
H the sum over k is dominated by the largest value of k at large enough r. The 

new element here is the sum over m; this summation is dominated by small m when 

4 2 

vT -C When a; <^ es, this condition automatically follows whenever vT -C 1. 
Consequently, naive perturbation theory is always good for times small compared to 
the inverse black hole temperature, in this regime. 

We emphasize that naive perturbation theory is always good at times of order 
St. Over such time scales (and for r ~ rj;/) we note that the sum over n and k are 
weighted by (this is as in section ^ while the sum over m is weighted by e^(-^)^. 
Note that the weighting factor for the sum over m is smaller than the weighting 

2 

factor for the sum over, for instance n, provided x <^ ea. It follows that our naive 
perturbation theory represents a weak departure from the black brane formation 
solution of section ^ when x -C es. 

2 2 

Now let us turn to the the parameter regime x ^ es. In this regime rj^R ~ 
so that black holes that are formed in the collapse process are always small in units 
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of the AdS radius. At times that are larger or of order 6t, the sum over m and n 

2 

are dominated by their smallest values provided "I ^ Making the replacement 
e = x'^e^ , this condition reduces to f <C ej-St which was exactly the condition for 
applicability of naive perturbation theory in flat space in section ^. The new element 
here is the sum over k. k is zero in section ^ and the sum over k here is dominated 
by = near r = th for <^ a condition that is automatically implied by 
-| -C Note, however, that, as in the previous paragraph, the sum over k is always 
dominated by the largest value of k at sufficiently large r. This reflects the fact that 
AdS space is never well approximated by a flat bubble at large r. Finally, specializing 
to V of order 5t and r ^ r^, the sum over n and m are each weighted by % ~ 4- 
while the sum over k is weighted by e^. In particular naive perturbation theory is 
good at times of order 6t provided x <^ y/e. 

Let us summarize in broad qualitative terms. Naive perturbation theory is a 

2 

good expansion to the true solution when vT <^ 1 for ^ <^ ^. In particular, this 
condition is always obeyed for times of order 6t when x <^ ^/e. 



4.5 Explicit results for naive perturbation theory 

As we have explained above, the functions 0, / and g may be expanded in an 
expansion in e as 



0(r, v) = e0i(r, v) + e^03(r, v) + 0{ef 

f{r, v) =rR{l + e'Mr, v) + e'U{r, v) + 0{eY) 

g(r, v) =r'^ + ^ + e^5(2(r, v) + eV(r, v) + 0(e)^ 



(4.17) 



Moreover the functions <^2n+i; fn and gn niay themselves each be expanded as a sum 
over two integer series (see ( [4.16|) ). The sum over k runs over a flnite number of 
values in ( [4.16| ) and we will deal with this summation exactly below. However the 
sum over m runs over all integers, and is computatble only after truncation to some 
flnite order. This truncation is justifled as the sum over m is effectively weighted 
by a small parameter as explained in the section above. In this section we present 
exact expressions for the functions 0i, g2 and /2, and expressions for 03, and (74 to 
the flrst two orders in the expansion over the integer m (this summation is formally 
weighted by 4); 
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The solutions are given as 



f2{r, v) 
92{r,v) 



12 



Wo C2{V) 



4 

K{v) 



1 



+ 



r_^{3K{t)-Mtr)dt ^ j:oodtii-oodt,(3K{t,)-Mt2r) 



12r4 



12r3 



h{r, v) 
94{r,v) 



\ 384r4 32r4 



V 96r4 120r5 >/ 



1 



^i(^) 15^3(^^) + '^A2(v) - 

24^2 48^3 



240r^ 



(4.18) 



Where 



dt 00 ( -C'2(t) + 0000 



Aiiv) = 0o(^^) / dh [' dt2 (-3K{t2) + 4it2)) 

J —00 J —00 

A2{v) = 0o(^) / dt (-3K{t) + ^l{t) 



(4.19) 



C2{V) = -\ 

1 



'o(t) 



dt (poit) 



■00 

V 



+ 0oW 



(4.20) 



4.6 The solution at late times 



t2 



dti 0o(ti) / dt2 / dt^ 3K{t3) - Mhr + C> 



R 



As in previous sections, our solution is normalizable (unforced) for v > St. Naive 
perturbation theory reliably establishes the initial conditions for this unforced evo- 
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lution at V = 6t. To leading order, this evolution is given by global AdS black hole 
metric with M = C2i6t) (see ( pTD ), perturbed by (f){6t) = see (|]T8]). As in the 
previous two subsections, the qualitatively important point is that this represents a 
small perturbation about the black hole background. Moreover, it follows on general 
grounds that perturbations in a black hole background in AdS space never grow un- 
boundedly (in fact they decay) with time. Consequently, we may rehably conclude 
that our spacetime takes the Vaidya form ( 4.1(J| ) at all times to leading order in the 
amplitude expansion. 

In order to determine an explicit expression for the subsequent dilaton evolution, 
one needs to solve for the linear, minimally coupled, evolution of a ^ initial condition 
in the background of global AdS with a Schwarzschild black hole of arbitrary mass. 
As in the previous two sections, the linear differential equation one needs to solve 
appears to be analytically intractable, but could easily be solved numerically. We 
will not, however, attempt this evaluation in this paper. 



5. Discussion 



In this paper we have used the AdS /C FT correspondence to determine the response 
of a conformal field theory, initially in its vacuum, to a low amplitude perturbation 
by a source coupled to a marginal operator. When the CFT in question lives on 
responds to the perturbation by the source by thermalizing into a plasma 
type phase. On the other hand, when the CFT in question lives on a sphere it 
either thermalizes into a plasma type phase or settles down into a glueball type 
phase depending on the details of the perturbation procedure. In this paper we have 
demonstrated that, to leading order in the amplitude expansion, the dual description 
of the thermalization into a plasma type phase is a spacetime of the Vaidya form. 
In odd boundary field theory dimensions the Vaidya metric reduces exactly to the 
uniform black brane metric in the causal future v > 6t at the boundary. As was 
discussed in detail in section plj], for many purposes our system behaves as if it has 
thermalized instantaneously. 

In this paper we have only studied solutions with a high degree of symmetry; for 
instance, solutions that maintain spatial translational invariance. The solutions of 
this paper may prove to be a useful starting point in describing the response of the 
field theory to a forcing function that breaks this symmetry, provided the scale of 
spatial variation of the forcing function is large compared to the inverse temperature 
of the black brane that is set up in our solutions. Consider, for example, the Einstein 
dilaton system studied in section g perturbed by the non normalizable part of a small 
amplitude dilaton field that takes the form (j)o{v,x). Let us further assume that the 
length scales for spatial variation L(x) in 0o are all large compared to (the 
inverse of the temperature of the black brane that is eventually formed). As e ^ 1 
this implies, in particular, that L{x) ^ 6t{x). We expect the resultant thermalization 
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process to be described by a dual metric of the form 



ds" = 2drdv - ( - ^^^^ ] dv' + rHx]. (5.1) 



where 



M{v,x) =C2{v,x) + 0{e^) 

1 r 

C2iv,x) = --J dt(j)o{t,x)(j)o{t,x) ^^2) 



oo 

oo 



C2{x) = \ I dt(f)o(t,x)(f)o{t,x) 

i.e. to be approximated tubewise by the solutions described in this paper. The 
metric ( |5.1| ) will then be corrected in a power series expansion in two variables; e 
(as in this paper) and a spatial derivative expansion weig hted by The last 

expansion should reduce to the fluid dynamical expansion at late times. Indeed, at 
t = St, the metric described in ( |5.1| ) is dual to a locally thermalized conformal fluid, 
everywhere at rest, but with a space varying energy density C2{x). The evolution 
of this fluid after v = 6t will simply be governed by the Navier Stokes equations 
of fluid dynamics; the metric dual to the corresponding flow was determined in 
I, I, I, I, I, I, 0, i I, H, I, O, H, ll, [H, m 0, |18|, [H. Provided we can 
solve the relevant fluid dynamical equations, we have a complete description of the 
evolution of our spacetime for all v. 

The gravitational solutions presented in this paper appear to be qualitatively 
different, in several ways, for odd and even d (see Appendix |B.ip . This suggests 
that the equilibration process at strong coupling is qualitatively different in odd and 
even dimensional field theories. At leading order in amplitudes, equilibration takes 
place faster in in field theories with odd spacetime dimensions as compared to their 
even dimensional counterpart. It would be interesting to find a direct field theory 
explanation of this fact. 

In this paper we have investigated the response of an AdS space to a marginal, 
non normalizable deformation of small amplitude. It would be natural to extend our 
work to determine the response of the same space to a relevant or irrelevant non 
normalizable marginal deformation of small amplitude. More ambitiously, one could 
also hope to explore the response of the system to large amplitude deformations, 
perhaps using a combination of analytic and numerical techniques (see |]S3|). We 
leave these issues to the future. 

In this paper we have constructed several solutions to bulk equations in a per- 
turbative expansion. It is natural to wonder whether the dynamical processes we 
have constructed in this paper are stable to small fluctuations, when embedded into 
familiar examples of the AdS /C FT correspondence. We will not address this ques- 
tion in detail in this paper; in these paragraphs we simply address the question of 
when the end point of the collapse processes, studied in this paper, are stable. 
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In this paper we have described time evolutions that end up in big black holes, 
small black holes (big and small compared to the AdS radius) and a thermal gas 
in AdS. To the best of our knowledge, large uncharged AdS black holes are stable 
solutions in every familiar example of the AdS/ C FT correspondence. The AdS 
thermal gas has a potential instability, the Jeans instability, which is triggered at 
energies at or larger than a critical density of order unity |7^, in the units of our 



paper. However the collapse situations described in this paper, that end up in a 

2 d-2 

thermal gas, do so at energies of order ^ -C e^-'^ ^ 1. We conclude that the 
thermal gases produced as the end point of collapse, in our paper, are also stable. 

Small black holes in AdSd+i x X, on the other hand, are usually unstable to a 
Gregory Lafiamme type clumping on the internal manifold X ||8^ . Consequently, the 



evolutions leading to small AdS black holes, constructed in this paper, are necessarily 
unstable when embedded into familiar examples of the AdS /C FT correspondence^^. 
Note, however, that the time scale associated with this Gregory Lafiamme instability 
is R^th where th is the Schwarszchild radius of the small black hole. Assuming 

2 

Rvh <^ 1 (so that the black hole that is formed is really small), ^^^^ ~ (^^T~") ^• 
It follows that, in the limit x'^^^ ^ e, studied in this paper, the black hole formation 
processes discussed in this paper occur over a time scale much smaller than that 
of Gregory Lafiamme instability. In other words the thermalization to small black 
holes (in the perturbative regime described in this paper) is described by a two 
stage process. In the first stage the solution is well described by the Vaidya metric 
(plus corrections) described in this paper. The second stage describes the evolution 
of an almost completely formed small AdS black hole perturbed by the Gregory 
Lafiamme instability. This black hole will then undergo the Gregory Lafiamme type 
transition in the usual manner. In other words the perturbative solutions presented 
in this paper correctly describe the process of small black hole formation even when 
embedded in AdS /C FT type situations in which this black hole is unstable. 

As we have explained above, a CFT on the sphere can respond to a forcing func- 
tion by settling down into either of its two available phases. It appears that the space 
of possible forcings is divided into two regions, by a critical surface of unit codimen- 
sion. On either side of this critical surface, the forcing drives the system into different 
phases. This critical surface occurs at x'^~^ ~ e^. Ignoring the issue of the Gregory 
Lafiamme instability for a moment, the local (small r) form of this gravitational solu- 
tion precisely on this critical surface is simply Choptuik's critical solution in i?'^"^'^, 
which is known to display surprisingly robust universal behavior characterized by a 
universal, nakedly singular solution and universal critical exponents. Note, however, 
that the Gregory Lafiamme instability generically cannot be ignored in the neigh- 
bourhood of the critical surface. In the neighbourhood of this transition black holes 
that are formed near criticality have arbitrarily small Schwarzscild radius, and so 
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trigger the Gregory Laflamme instability over very short time scales. Consequently, 
in order to access the universal i?*^"^'^ Choptuik behavior in field theory one would 
have to tune initial data with exponential accuracy. 

It is of course true on general grounds that the solutions described in this paper 
fall into two classes distinguished by an order parameter ( the presence of a hori- 
zon at late times). An important question about the transition between these two 
behaviours is whether it is continuous or discontinuous over time scales small com- 
pared to In situations in which Gregory Lafiamme instability occurs there is a 
very special submanifold of this transition manifold; the submanifold on which Gre- 
gory Lafiamme instabilities are precisely tuned away. On this submanifold we know 
that the relevant gravitational solutions are those of Choptuik collapse in and 
so are continuous (second order) and singular. We do not konw if this singular sec- 
ond order behaviour persists away from this special point. The investigation of these 
issues, as well as the study of the smoothing out of singularities on this manifold by 
finite fiuctuations, is potentially interesting area of future research. 
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6. Appendices 

A. Trans lationally invariant graviton collapse 

In sections |^ and ^ above we have studied the collapse triggered by a minimally 
coupled scalar wave in an asymptotically AdS background. Our study was, in large 
part, motivated by potential applications to CFT dynamics via the AdS /C FT corre- 
spondence. From this point of view the starting point of our analyses in e.g. section 
1^ has a drawback as not every bulk system that arises in the study of the AdS /CFT 
correspondence, admits a consistent truncation to the theory of gravity coupled to a 
minimally coupled massless scalar field. 
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On the other hand, every two derivative theory of gravity that admits AdS space 
as a solution admits a consistent truncation to Einstein gravity with a negative 
cosmological constant. Consequently, any results that may be derived using the 
graviton instead of dilaton waves, applies universally to all examples of the AdS / CFT 
correspondence with two derivative gravity duals. In this section we study a situation 
very analogous to the set up of section 0, with, however, a transverse graviton playing 
the place of the dilaton field of section 0. All the calculations of this section apply 
universally to any CFT that admits a two derivative gravitational dual. 

While the equations that describe the propagation of gravity waves are more 
complicated in detail than those that describe the propagation of a massless min- 
imally coupled scalar field, it turns out that the final results of the calculations 
presented in this subsection are extremely similar to those of section |^. We take 
this to suggest that all the qualitative results of sections | and ^ would continue to 
qualitatively apply to the most general approximately translationally invariant grav- 
itational perturbations of Poincare Patch AdS space or approximately spherically 
symmetric gravitational perturbation of global AdS space. If this guess is correct, 
it suggests that the qualitative lessons learnt in this paper have a wide degree of 
applicability. 

In this section we restrict our attention to the simplest dimension d = 3. It 
should we possible, with some additional effort, to extend our results at least to all 
odd d, and also to work out the corresponding results for even d. We leave this 
extension to future work. 

The set up of this Appendix is very closely analogous to that employed by Yaffe 
and Chesler in [^. The main differences are as follows. Yaffe and Chesler worked 



in d = 4; they numerically studied the effect of a specific large amplitude non nor- 
malizable deformation on the gravitational bulk. We work in c/ = 3, and analytically 
study the the effect of the arbitrary small amplitude deformation on the gravitational 
bulk. 

A.l The set up and summary of results 

In this section we study solutions to pure Einstein gravity with a negative cosmolog- 
ical constant. We study solutions that preserve an R"^ x Z2 x Z2 symmetry. Here 
denotes the symmetry of translations in spatial field theory directions, while the two 
Z2S respectively denote the spatial parity flip and the discrete exchange symmetry 
between the two Cartesian spatial boundary coordinates x and y. 

As in section ^, our symmetry requirements determine our metric up to three 
unknown functions of v and r. With the same choice of gauge as in section |^, our 
metric takes the form 

ds"^ = -2 dv dr + g{r,v) dv^ + f{r,v){dx^ + dy"^) + 2r'^h{r,v)dx dy (A.l) 
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The boundary conditions on all fields are given by ( p.9|) under the replacement 
4>{r,v) — > h{r,v) and (po{v) ho{v). Here ho{v) gives the boundary conditions 
on the off diagonal mode, g^y, of the boundary metric. hQ{v) is taken to be of order 
e. Physically, our boundary conditions set up a graviton wave, with polarization 
parallel to the spatial directions of the brane. 

As in section 0, in order to solve Einstein's equations with the symmetries above, 
it turns out to be sufficient to solve the three equations E^, Eq and E^y (see ( |2.10| )) 
(plus the energy conservation condition rEec at one r). 

As in section ^ it is possible to solve these equations order by order in e. We 
present our solution later in this section. To end this subsection, we list the principal 
qualitative results of this section. We are able to show that 

• The boundary conditions described above result in black brane formation for 
an arbitrary (small amplitude) source function hQ{v). 

• Outside the event horizon of our spacetime, we find an explicit analytic form 
for the metric as a function of ho{v). Our metric is accurate at leading order 
in the e expansion, and takes the Vaidya form ( p..l|) with a mass function 

M{v) = -\ j dthoho (A.2) 

• In particular, we find that the energy density of resultant black brane is given 
by 

1 f°° 

M ^ -E2 = - dthl (A.3) 

2 J -00 

2 

Note that E2 ~ 

• As this leading order metric is of the same form as that in the previous sub- 
section, the analysis of the event horizons presented above continues to apply. 
In particular it follows that singularities formed in the process of black brane 
formation are always shielded by a regular event horizon at small e. 

• Going beyond leading order, perturbation theory in the amplitude e yields sys- 
tematic corrections to this metric at higher orders in e. We unravel the structure 
of this perturbation expansion in detail and work out this perturbation theory 
explicitly to fifth order at small times. 

A.2 The energy conservation equation 

As we have explained above, the equations of motion for our system include the 
energy conservation relation, in addition to the one dynamical and two constraint 
equations. The form of the dynamical and constraint equations is easily determined 
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using Mathematica-6; these equations turn out to be rather lengthy and we do not 
present them here. In this section we content ourselves with presenting an explicit 
form for the energy conservation equation. As in section ^, it is possible to solve for 
the functions -, and h in a power series in -. This solution is simply the Graham 
Fefferman expansion. To order ^ (relative to the leading result) we find 



/ 



f{r,v) = r 
g{r, v) = 

/ 



ho 



1 + 



+ 



+ 



V 



V 



+ 



1 








{l + 3hl) 


ho 


' -Ahoi-l + hl)d'^ho 



M{v) 



+ 



J 



h{r, v) 



hohi 



, ,ho 4(-mg) a{v) 

fto H 1 o 1 5 r 



v 



(A.4) 



where the parameters M and a are constrained by the energy conservation equation 
^° ^ [ + 3M{v)ho (-1 + hlY-3 (-1 + hlYa 

(A.5) 



M 



4 (-1 + hi) h„h„alh„ + (-1 + A^'aj/io + (1 + ihi) 





3 " 


ho 





24 



In the perturbative solution we list below, we will find that a ~ 0{t^). It follows 
that, to order C(e^), the function M{y) is given by ( |A.2| ). 

A. 3 Structure of the amplitude expansion 

As in subsection ^ we set up a naive amplitude expansion by formally replacing ho 



24 



The stress tensor is given by 
Ttt = M 

-Lxx — -Lyy — ^ 



T — — 



2{-l + hl) 

+ {-i+hiyhi + {i+3hi) 



- 3 (-1 + K) ai^) - 4 (-1 + K) hohodtho 

3' 



(A.6) 



Using these relations, it may be verified that ( [A.5| ) is simply a statement of the conservation of the 
stress tensor. 
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with eho and then solving our equations in a power series in e. We expand 

oo 

f{r,v) = J2^yn{r,v) 

n=0 

oo 

g{r,v) = J2^''9n{r,v) (A.7) 

n=0 

CX) 



n=0 

with 



/o(r,t;) = r, 5(o(r,f) = r^ ho{r,v) = 0. (Ai 



The formal structure of this expansion is identical to that described in section ^.5[ 
in particular /„ and Qn are nonzero only for even n while hn is nonzero only for odd 
n. At first order we find 

h,(r,v) = ho{r,v) + ^^^^^ (A.9) 

r 

which then leads to simple expressions (see below) for /2 and g2- In particular hi 
and /2 vanish for v > 6t while (72 = M/r for > 6t. 

Turning to higher orders in the perturbative expansion, it is possible to induc- 
tively demonstrate that for n > 1 

• 1. The functions hn, g-n and /„ have the following analytic structure in the 
variable r 



h2n+i{r,v) = 2^ 

k=2 
2n+2 N 

/2.+2(r,^)=r^^ (A.IO) 

k=2 



n 



92n+2{r,v) = 



k=l 



• 2. The functions /?-2n+i('^)! /2n+2('^) ^'^d 5'2n+2(^) each functional of /io('i^) 
that scale like A~'^ under the scaling v —>■ Xv. 

• 3. For V > 6t these functions are all polynomials in f of a degree that grows 
with n. For example, the degree of /i2n+i most 3n — k. 

As in the section ^, this structure ensures that naive perturbation theory is good 
for times v <^ M^, but fails for later times. As in section (^, the correct perturbative 
expansion uses the Vaidya metric ( p,.l| ) as the zero order solution. 
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A. 4 Explicit results up to 5th order 

At leading order we have 



hi{r,v) = ho{v) + 



ho 



f2ir,v) 
92{r,v) 



hn 



E2{v) ^ 1 
r 4 



+ hodlhQ 



(All) 



At next order 
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he 



h{r,v) 



hl{v) 



he 
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D{v)ho{v) ho 
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128r3 
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/i4(r, = 
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■oo 



(A12) 



Finally at the next order 



h^{r,v) 
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24r3 
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(A.13) 
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where 

D2{x) 
D^{x) 

DJx) 



-hoixf[d.Mf 

Ei{x)dxho + '^D{x)ho{x)dxhQ + E2{x)hQ{xY d^ho 
\ 5D{x) [d.hof + 15E2{x)ho{x) [d.hof + 15D{x)ho{x)dlho 
^ 18D{x)E2{x) + 5E2{x) [dMf + 7D{x)ho{x)dlho 



(A.14) 



and (this follows from energy conservation) 
1 • 



E2 = ^hod% 



3E2iv)hoiv) 



ho 



+ Aho{v)hodtho + 2hldlh 



(A.15) 



It follows in particular that the the 'initial' condition for normalizable evolution 
ai V = 5t is given, to leading order, by 



h{r, St) 



—00 \j —00 



dy {dyhodyho) dxho{x)d 



X 



(A.16) 



This initial condition is of order 



i.e. of order 4t where r 



This demon- 
strates that, for V > 6t, our solution is a small perturbation about the black brane 
of energy density E2. 

A. 5 Late Times Resummed perturbation theory 

To leading order, the initial condition for the normalizable evolution of resummed 
perturbation theory for the field h{r, v) is given by 



HSt) = ^ E2{x)dxho dx^ = 



him 



Now, at the linearized level the equation of motion for the function h is simply the 
minimally coupled scalar equation. It follows that the subsequent evolution of the 
field h is simply given by 



M Mi 



(A. 17) 



where the universal function if) was defined in section ^. As in section H, this per- 
turbation is small initially, and at all subsequent times, justifying the resummed 
perturbation procedure. 
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B. Generalization to Arbitrary Dimension 

B.l Translationally Invariant Scalar Collapse in Arbitrary Dimension 

In this subsection we will investigate how the results of section ^, which were worked 
out for the special case d = 3, generalize to c? > 3. The mathematical problem we 
will investigate in this Appendix was already set up in general d in subsection |2.1| . It 
turns out that the dynamical details of collapse processes in odd and even dynamics 
are substantially different, so we will deal with those two cases separately. 



B.1.1 Odd d 

The general structure of the solutions that describe collapse in odd > 5 is similar in 
many ways to the solution reported in section ^ The energy conservation equations 
may be studied via a large r Graham Fefferman expansion closely analogous to that 
described in section 0. The functions / and g may be expanded at large r as 

n=0 

\n=0 / 

g{r,v)=r [}_^^:^ 

\n=0 / 

For n < d ~ 1 the equations of motion locally determine A^lv), A^{v) and Ag{v) 
in terms of (poiv)- Each of these functions is a local expression (of n^^ order in v 
derivatives) of (f)o{v). However local analysis does not determine Ag[v) = M{v) and 
A'^{v) = L{v) in terms of (f)o{v)- M{v) and L{v) are however constrained to obey an 
energy conservation equation that takes the form 

M = k(j)L{v) + local (B.2) 



where is a constant and 'local' represents the a set of terms built out of products of 
derivatives of (f)o{v) that we will return to below. As in c? = 3, L{v) = 0{e^), so the 
first term in (p.2| ) does not contribute at lowest order of the amplitude expansion of 
interest to this paper. The local terms in this equation ( p.2|) are easily worked out at 
lowest order, O(e^), in the amplitude expansion, and we find M(t>) = C2(f) + 0{e^) 
with 



-)(i-2 



C2iv) 



(¥)' 




d+3 



d-3 



(B.3) 
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od-l / (<t±)\ \ ^ roc 2 2 



the generalization of (|2.2CI|) and ( p.21| ) to arbitrary odd d. ( p.4|) gives the leading 



order expression for the mass of the black brane that is eventually formed at the end 
of the thermalization process. 

Let us now turn to the naive amplitude expansion in arbitrary odd d. The first 
term in this expansion, 0i is easily determined and we find 



Equations ( p.l2| ) then immediately determine /2 and g2. Turning to higher orders, 



it is possible to demonstrate that 

• 1. The functions 4>2n+i, fl'2n and /2„ have the following analytic structure in the 
variable r 




(2n + l)(d-l) 

02n+l(^) 



2 

(2n4 
T 2 



(2n + l)(d-l) 



■k 



j,n{d—l)—k 

k=0 



(B.6) 



^ / N / n{d-l)-g(n) ^ , , 

92n{r,V) - + 1^ ^n{d-l)-k 



\ k=0 

where 

p{n) = d, {2n + l>d), p{n) = 2n + I {2n + l<d), 

fin) = d, {2n > d), f{n) = 2n {2n < d), 
g{n)=d-l, (2n>rf-l), ^(n) = 2n - 1 {2n < d). 

2. The functions 02n+i(^)? /Inl^) 92ni'^) each functionals of 0o(f ) that 
scale like A~'^ under the scaling v Xv. 

3. For V > 6t f2 = U = 0, g2 = and = Further, effectively, 
p{n) = d, f{n) = 2d and g{n) = 2d — 1 for v > 6t (all additional terms present 
in ( p.6| ) vanish at these late times). Moreover the functions 02n+i(^)' /2n(^) 
and 5'2n('^) are all polynomials in v whose degrees are bounded from above by 
n + k — l,n + k — 3 and n + k — 4 respectively. 
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As in (i = 3, the polynomial growth in v of the coefficients of the naive per- 
turbative expansion invalidates this expansion for large enough v. More specifically, 

2 

the sums over k and n in the expressions above are weighted by rv and -^r respec- 

2 

tively. In the neighborhood of the horizon, r ^ th ^ T ^ each of these sums is 
effectively weighted by the factor vT. Consequently, naive perturbation theory fails 
at times large compared to the inverse temperature of the brane. At times of order 

2 

8t and for r th the sum over k and n are each weighted effectively by e^. More 

2 

generally, naive perturbation theory is good at times of order 5t provided r5t ^ e<*-i , 
a condition that is satisfied at the event horizon. 

As in (i = 3 the IR divergence of the naive perturbation expansion has a sim- 
ple explanation. Even within the validity of the naive perturbation expansion, the 
spacetime is not well approximated by empty AdS space, but rather by the Vaidya 
metric ( |1 . 1| ) . The naive expansion, which may be carried out with comparative ease 
up to f = may be used to supply initial conditions for the subsequent unforced 
normalizable evolution for resummed perturbation theory. For v > 6t, the spacetime 
metric is given, to leading order, by the Vaidya form ( [L.lD , with C2{v) given by the 
constant C2 listed in (p.4| ) 

Consequently, the spacetime metric for v > St is the black brane metric with 

2 

temperature of order perturbed by a propagating field and consequent space- 
time ripples. The initial conditions at v = 6t, that determine these perturbations at 
later times, are given to leading order in e (read off from the most small r singular 
term in 03 ) as 



0(r,f) 
where 



A 



3(d-l) 
f 2 



2{d-2) J_ 



00 



(¥)' 
{d-l)\ 



3 

d-l \ 2 / d+1 



- 1 



In terms of the normalized variable x = and y = vM d this initial condition 

M d 

takes the form 

3 

m ~ (B.8) 

X 2 

It follows that the solution a.t v > 6t is (in the appropriate x, y coordinates) 

3 

an order perturbation about the uniform black brane. The coefficient of this 
perturbation is bounded for all y, and decays exponentially for large y over a time 
scale of order unity in that variable. The explicit form of the solution for 0, for 
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V > 6t, may be obtained in terms of a universal function, 4^d{x, y) as in section ^ 
The equation that we need to solve is 




(B.9) 



Figure 5: Numerical solution for the dilaton at late time in d = 5 

As in section 0, this universal function appears to be difficult to obtain analyt- 
ically, but is easily evaluated numerically. As an example in Figure ^ we present a 
numerical plot of this function in (i = 5. As in section ^ we find it convenient to 
display the numerical output for the function ip5{\,y) over the full exterior of the 
event horizon, w G (0, 1). 

In figure ^ we present a graph oiifj^^-^^y) (i.e. as a function of time at a fixed 
radial location) Notice that this graph decays, roughly exponentially for v > 0.5 
and that this exponential decay is dressed with a sinusodial osciallation, as expected 
for quasinormal type behavior. A very very rough estimate of this decay constant 
is provided by the equation ui using the equation 'f'K^'l.l = e'^-^"^^ which gives 

''/'s (, qty I ■ j 

^^In order to obtain this plot, as in||, we worked with the redefined field XsC^: v) = (l~")V'5(:^j 2/) 
and imposed Dirichlet boundary conditions on this field at u = and u = 0.999999. We also imposed 
the initial conditions X5 = (0.999999— The Figure |5| was outputted by Mathematica-6's partial 
differential equation solver, with a step size of 0.0005 and an accuracy goal of 0.001. 
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uji ^ 8.2T (here T is the temperature of our black brane given by T = This 
number is the same ballpark as the decay constants for the first quasi normal mode 
of the uniform black brane reported in (unfortunately those authors have not 
reported the precise numerical value for c? = 5) . 



B.1.2 Even d 



In our analyses above we have so far focused attention on odd d (recall that d is 
the spacetime dimension of the dual field theory). In this subsection we will study 
how our results generalize to even d. While all the broad qualitative conclusions of 
the odd d analysis plausibly continue to apply, several intermediate details are quite 
different. The analysis of all equations is more difficult in even than in odd dimensions. 
In this appendix we aim only to initiate a serious analysis of these equations, and 
to carry this analysis far enough to have a plausible guess for the behavior of our 
system. We leave a systematic analysis of these equations to future work. 

The qualitative differences between even and odd d show themselves already in 
the Graham Fefferman expansion. We illustrate this by working out this expansion 
m d = A. In this dimension the expansion of /, g, at large r take the form 
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f( \= _(0o£_M :-3(0o)' + 20000 -aVo)' 
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-190o(0o)3 - 144OL(^)0o - 18(9,^0^,00 + 450o0o 
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(B.IO) 



12'" ' r2 24r2 

■y^ 



log(r) a>o - 2(0o)'920o 



16r4 

The energy conservation equation is 

M = (4O0o(0o)' - 192L(t;)0o - 179>o0o + 60o0o) (B.ll) 
and at quadratic order in e we have 
M{v) = C2{v) + O(e^) 

1 r / o ••• \ (B.12) 

C2(t^) = — y rft(-192L(t)0o-17 0o0o + 6920o0o) 

Unhke in even dimensions, it turns out that in odd dimensions L{v) is nonzero 
at order e. This is fortunate, as all the local terms in (|B.12| ) are total derivatives, 
and so vanish when v is taken to be larger than 5t. The full contribution to the mass 
of the black brane that is eventually formed from our collapse process arises from 
the term in ( p.l2|) that is proportional to L{v). As a consequence, the mass of the 
eventual black brane is not determined simply by Graham Fefferman analysis, but 
requires the details of the full dynamical process. These details may be worked out 
at lowest order in the e expansion, (see below) and we will find 

L{v) = d% + ^ £ dthgiv - t) dUoit) + 0{e') (B.13) 

Plugging into (p.l2| ) we find that C2(f) reduces to the constant C2 for v > 6t, and 
we have 

1 f'^ 

C2 = -— dhdt2 {dlMh) log(ti - t2)e(ti - t2)dlMh)) (B.14) 
-i-^ J-00 

Let us now turn to the amplitude expansion of our solutions. We will work this 
expansion out only at leading order; already the leading order solution turns out to 
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have qualitative differences (and to be much harder to determine and manipulate) 
than the corresponding solution in odd d. 

Recall that (pi ( |B.5| ) is extremely simple when d was odd. To start with, the 
solution is local in time, i.e. 0i(r, wq) is completely determined by the value, and a 
finite number of derivatives, of 0o(wo)- Relatedly (/>(r, f) has a very simple analytic 
expression in r; it is a polynomial in ^ of degree In even d, on the other hand 
the dependence of 0i(r, on (poiv) is not local in time. Relatedly, the expansion of 
0i(r, f) in a power series in ^ has terms of every order in K Explicitly we find 




(B.15) 



where cosh 9 = 1 + x 



Note that the function h[x) admits the following large x expansion 
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^ {d-k)k\{d- 



\m=l 



+ 
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^ln(2x) + 0(l!i£ 



2''(i!) 
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B.16) 



The fact that h{x) grows (rather than decays) with x may cause the reader to worry 
that (f)(r, v) blows up at large v. That this is not the case may be seen by noting 
that v'^ d'^'^^ (po may be rewritten as a sum of total derivatives when k < d + 1 and so 
integrates to zero when v > 6t (in general it integrates to a simple local expression 
even for v < 6t). Explicitly, plugging ( |B.16|) into ( B.15 ) and integrating by parts we 
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find that 0i(r, v) has the following large rt behavior 

d 



Ai{v) B{v)\n{r) 



+ 0{ 



Inr , 
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(B.17) 



+ 0( 
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■ ln(2rt) 



d+l 
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(where the functions Ai{v) and B{v) are defined by this equation). On the other 
hand at small x we have 



h{x) 



{2xy 
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(B.18) 



from which it follows that 

, , . 1 1 



)ilr,f I 



dt{2{v - t))^dt'Mt) + Oi^), (B.19) 



T 2 



an expression that is valid at small rv. Note, in particular, that for 6t <^ v, (p.l9| ) 
reduces to 



1 



Mr,v) = + 0(4.) + 



(B.20) 



In particular this formula determines the behavior of the field 0i in the neighborhood 
of the event horizon r// ~ T for times that are large compared to 6t but small 
compared to T^^. 

The functions /2 and g2 are easily expressed in terms of the function (f)Q. We find 



f2ir,v) 
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r / {dpcf>ifdp 



p\dp<P,fdp 



2{d-l) 

g2{r, v) = - {2dj2{r, v) + {d ~ 2)rh{r, v) + r^ajslr, v)) (B.21) 
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„d-2 



- 64 - 



The function D2{v) is determined by the requirement that the coefficient of in 
the large r expansion of g2{r,v) is —C2{v) (see ( P.12|) ); in particular, for v > 6t, 
D2{v) = C2{v). At small r and for v > 6t 

' ' 2{d - l)id - 2){d - ?,)rd-^ ^^rd-^> 



, X C2 d^K^iv) 1 



r^-^ {d - l){d - 2){d - 3)r<^-^ V'^-S' 



(B.22) 

{d + l){d-l)\ J_ 

2^?J>^ 1 

i.t^ (d+ind-iy. ^"■^"■> 

We would like to draw attention to several aspects of these results. First note 
that 0i(r, t>) is small provided (r5t)~ ^ e. Consequently, we expect a perturbative 
analysis to correctly capture the dynamics of our situation over this range of coor- 
dinates; note that this is exactly the same estimate as for odd d. Next note that 
the maximal singularity, at small r, in the functions /2 and g2, are both of order 
this is the same as the maximal singularity in the analogous functions in odd 
d (see the previous subsection). As the function g^ir^v) = r^, it follows, as in the 
previous function, that our spacetime metric is not uniformly well approximated by 
the empty AdS space over the full range of validity of perturbation theory. Over this 
entire range, however, it is well approximated by a Vaidya type metric, where the 
mass function for this metric is given at leading order by the coefficient of —-^z^ in 
g2{r, v) above. 

Unlike the situation in odd dimensions, the leading order mass function M(f), in 
the effective Vaidya metric, is not given simply by C2{v). In particular, when v ^ 5t 
we have from ( p.22|) that 



C2 - M{v) f5t\ 



C2 \v 



d+2 



In other words, the leading order metric for the thermalization process, in even d, 
is not given precisely by the metric of the uniform black brane for v > 6t. However 
it decays, in a power law fashion, to the black brane metric at times larger than St. 
As a consequence at times 5t <^v <^ the leading order metric that captures the 
thermalization process is arbitrarily well approximated by the metric of a uniform 
black brane. It follows that, while the spacetime described in this subsection does 
not capture the dual of instantaneous field theory thermalization (as was the case in 
odd d), it yields the dual of a thermalization process that occurs over the time scale 
of the forcing function rather than the much longer linear response time scale of the 
inverse temperature. 
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We will not, in this paper, continue the perturbative expansion to higher orders 
in e. We suspect, however, that the computation of 03 when carried through will 

3 

yield a term proportional to ^^a-i) that is constant in time. This term will dominate 

r — 2 — 

the decaying tail of 0i(r, f) at a time intermediate between St and and will set 
the initial condition for the late time decay of the field (over time scale T~^) as 
was the case in odd dimensions. It would be very interesting to verify or correct this 
guess. 

B.2 Spherically Symmetric fiat space collapse in arbitrary dimension 
B.2.1 Odd d 

The discussion of section |^ also extends to the study of spherically symmetric collapse 
in a space that is asymptotically flat R^'^ for arbitrary odd d. In this section we will 
very briefly explain how this works, focussing on the limit y = ^ ^ 1. 
To lowest order in the amplitude expansion we find 

Mr,v)-l^2^ m! (^-m)! ^^'^^^ 

Here i^^v) is a function of time that we take, as usual, to vanish outside v G (0, 6t), 
and be of order ef{6t)~ , where e/ is a dimensionless number such that e/ ^ 1. As 
in section ^ the parameter that will justify the amplitude expansion will be ^. 

( [B.23| ) together with constraint equations immediately yields an expression for 



the functions /2 and g2. In particular, the leading large r approximation to g2 is 
given by 



M{v) 

M{v) = -j-j I dt 



(d-3) \ / jd+l) \ (7 _ 3 / (d-1) 



(B.24) 



Note that 0i ^ 1 whenever r~ <^ {6t)^^ef so we expect the amplitude ex- 
pansion to reliably describe dynamics over this range of parameters. As in section 
1^, however, g2 cannot be ignored in comparison to gQ = 1 throughout this parameter 
regime. As in section ^, this implies that our spacetime is well approximated by a 
Vaidya type metric rather than empty fiat space even at arbitrarily small ^. The 
mass function of this Vaidya metric is given by M{y) in (p.24| ). 

As in section |^ one may ignore this complication at early times v ^ th over 
which the solution is well approximated by a naive perturbation expansion that uses 
empty flat space as its starting point. It is possible to demonstrate that this naive 
expansion has the following analytic structure in the variables r and v 
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1. The functions $2n+i, and G2n have the following analytic structure in 
the variable r 



m=0 ' 
m=0 

m=0 



• 2. The functions $2n+i(^)? -^2n(^) ^^"^ ^2n(^) functionals of ) that 

scale like Xm'(2n+l)i^ ^m-n(d~3) ^m-n(d-3)-l ^^^J^^ ^J^g ^J^^ SCallug V ^ Aw. 

M(f ) scales like A^"*^ under the same scaling. 

• 3. For V > St the $2n+i('^) polynomials in v of degree < n + m — 1; F^{v) 
and are polynomials in v of degree < n + m — 3 and n + m — 4 respectively. 

It follows that, say, 0(r, t>), is given by a double sum 



0(r,t.) = ^$2„+i(r,t;) = ^ 



{2n+l)i^+m' 

n n,m=0 

Now sums over m and n are controlled by the effective expansion parameters ~ ^ 
(for m) and ^.^^^J^a^-i ~ — iir ~ ^ (fo^ '^i recall that in the neighborhood of the 

horizon r'jf^ ~ ((5t)'^-2e^). 

As in section ^, it follows that the naive perturbation expansion breaks down 
for times v ^ th- However this expansion is valid everwhere outside the event 
horizon at times of order 5t, and so may be used to set the initial conditions for a 
resummed perturbation expansion that uses the Vaidya metric as its starting point. 
For V > 6t the mass function of the Vaidya metric reduces to a constant. At long 
times our solution is given by a small perturbation around a black hole of mass M. 
This perturbation is best analyzed in the coordinates x = — and y = — 

M'i-2 A/ d-2 

In these coordinates the leading order tail of 0, at long times, is given by motion 
about a black hole of unit Schwarzschild radius perturbed by the field with initial 
condition 

1 



0(x,O) 



3(d-l) 3{d-l) -2_ 



The smallness of this perturbation justifies linearized treatment of the subsequent 
dynamics. 
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B.2.2 Even d 



We will not, in this paper, attempt an analysis of the spherically symmetric collapse 
to form a black hole asymptotically W^'^ for even d. Here we simply note that the 
leading order large e/ solution for (j)i{v) may formally be expressed as 



dco |^g(a;)e-(--) j (B.26) 

for any function q{u}) where Hn{x) is the n*^ Hankel function of the first kind, i.e. 

Using this expansion, it is easily verified that 0i(r, f) reduces, at large r, to an in- 
coming wave that takes the form -^^rr- The evolution of this wave to small r is 
implicitly given by ( [B.26| ). It should be possible to mimic the analysis of subsub- 
section |B.1.2| to explicitly express 0i(r, f) as a spacetime dependent Kernel function 



convoluted against il){v). In analogy with subsection |B.1.2| it should also be possible 



to expand g2{j', v) about small r. It is tempting to guess that such an analysis would 
reveal that the leading singularity in g2{r,v) scales like so that the metric is 
well approximated by a spacetime of the Vaidya form. We leave the verification of 
these guesses to future work. 

B.3 Spherically symmetric asymptotically AdS collapse in arbitrary di- 
mension 

It should be straightforward to generalize the analysis of section ^ to arbitrary odd d, 
and perhaps also to arbitrary even d. We do not explicitly carry out this generaliza- 
tion in this paper. However it is a simple matter to infer the various scales that will 
appear in this generalization using the intuition and results of subsections |B.1| and 
B.2|, and the fact that the results of global spherically symmetric AdS collapse must 



reduce to Poincare patch collapse in one limit and fiat space collapse in another. We 
have reported these scales in the introductionto section ^. 
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